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Abstract 

In this paper we study geometries on the manifold of curves. 

We define a manifold M where objects c £ M are curves, which we 
parameterize as c : S 1 — > 1R™ (n > 2, S 1 is the circle). Given a curve c, we 
define the tangent space T C M of M at c including in it all deformations 
h : S 1 -> R" of c. 

We discuss Riemannian and Finsler metrics F(c, h) on this manifold 
M, and in particular the case of the geometric H° metric F(c, h) — 
J \h\ 2 ds of normal deformations h of c; we study the existence of min- 
imal geodesies of H under constraints; we moreover propose a conformal 



1 Introduction 

In this paper we study geometries on the manifold M of curves. This manifold 
contains curves c, which we parameterize as c : S 1 — ► IR" (S 1 is the circle). 
Given a curve c, we define the tangent space T C M of M at c including in it 
deformations h : S 1 — > IR™, so that an infinitesimal deformation of the curve c 
in direction ft, will yield the curve c(9) + eh(0). This manifold M is the Shape 
Space that is studied in this paper. 

We would like to define a Riemannian metric on the manifold M of curves: 
this means that, given two deformations h,k £ T C M, we want to define a scalar 
product (h, k) c , possibly dependent on c. The Riemannian metric would then 
entail a distance d(co,ci) between the curves in M, defined as the infimum of 
the length Len(7) of all smooth paths 7 : [0, 1] — > M connecting Co to ci. We 
call minimal geodesic a path providing the minimum of Len(7) in the class of 7 
with fixed endpoints. 

A number of methods have been proposed in Shape Analysis to define dis- 
tances between shapes, averages of shapes and optimal morphings between 
shapes; some of these approaches arc reviewed in section §2. At the same time, 
there has been much previous work in Shape Optimization, for example Image 
Segmentation via Active Contours, 3D Stereo Reconstruction via Deformable 
Surfaces; in these later methods, many authors have defined Energy Functional 
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on curves (or surfaces) and utilized the Calculus of Variations to derive curve 
evolutions to minimize the Energy Functionals; often referring to these evolu- 
tions as Gradient Flows. For example, the well known Geometric Heat Flow, 
popular for its smoothing effect on contours, is often referred as the gradient 
flow for length. 

The reference to these flows as gradient flows implies a certain Riemannian 
metric on the space of curves; but this fact has been largely overlooked. We 
call this metric H° henceforth. If one wishes to have a consistent view of the 
geometry of the space of curves in both Shape Optimization and Shape Analysis, 
then one should use the H° metric when computing distances, averages and 
morphs between shapes. 

In this paper we first introduce the metric H in §2.1; we immcdiatly re- 
mark that, surprisingly, it does not yield a well define metric structure, since 
the associated distance is identically zero*- 1 ). In §4 we analyse this metric; we 
show that the lower-semi-continuous relaxation of the associated energy func- 
tional is identically zero (see §3.1 and 4.14); but we prove in thm. 4.20 that, 
under additional constraints on the curvature of admissible curves, the metric 
H° admits minimum geodesies; we propose in §4.6 an example that justifies 
some of the hypotheses in 4.20. We can then define in §4.5 the Shape Space of 
curves with bounded curvature, where the metric H° entails a positive distance. 
These hypotheses on curvature, however, are not compatible with the classical 
definition of a Riemannian Geometry. 

More recently, a Riemannian metric was proposed in [MM] for the space M 
of curves (see §2.1.3 here); this metric may fix the above problems; but it would 
significantly alter the nature of gradient flows used thus far in various Shape 
Optimization problems (assuming that one wishes to make those gradient flows 
consistent with this new metric). In this metric, distances measured between 
curves are defined using first and second order derivatives of the curves (and 
therefore the resulting optimality conditions involve up to fourth order deriva- 
tives); as a consequence, flows designed to converge towards these optimality 
conditions are necessarily fourth order, thereby precluding the use of Level Set 
Methods which have become popular in the field of Computer Vision and Shape 
Optimization. 

We propose instead in §5 a class of conformal metrics that fix the above 
problems while minimally altering the earlier flows: in fact the new gradient 
flows will amount to a simple time reparametcrization of the earlier flows. In 
addition the conformal metrics that we propose have some nice numerical and 
computational properties: distances measured between curves are defined using 
only first order derivatives (and therefore the resulting optimality conditions 
involve only second order derivatives); as a consequence, flows designed to con- 
verge towards these optimality conditions are second order, thereby allowing 
the use of Level Set methods: we indeed show such an implementation and a 
numerical example in §5.3. We also proposed in [YM04] a differential operator 
that is well adapted to the problem at hand: we review it here as well, in §5.1.1. 



MThis striking fact was first described in [Mum] 
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1.1 Riemannian and Finsler geometries 

Let M be a smooth connected differentiable manifold. ^ For any c £ M, let 
T C M be the tangent space at c, that is the set of all vectors tangent to M at c; 
let TM be the tangent bundle, that is the bundle of all tangent spaces. 

Definition 1.1 Let X be a vector space; a norm \ ■ \ satisfies 

1. \v\ is positive homogenous, i.e. 

A|w| = |A«| VA>0 , 

2. 

\v + w\ < \v\ + \w\ 
(that, by (1), is equivalent to asking that \v\ be convex), and 
3. \v\ = only when v = 0. 
If the last condition is not satisfied, then | • | is a seminorm. 

Definition 1.2 We define a Finsler metric to be a function F : TM — > R + , 
such that 

• F is lower semicontinuous and locally bounded from above, and, 

• for all c £ M , v i— > F(c, v) is a norm on T C M . ( 3 ' 

We will also consider it to be symmetric, that is, F(c, —v) — F(c,v). 

We will write |u| c for F{c, v), or simply \v\ when the base point c can be 
easily inferred from the context. 

We define the length of any locally Lipschitz path 1 - 4 ) 7 : [0, 1] — > M as 



Len7 = 




and the energy as 

E( 7 )= f \i(v)\ 2 l(v) dv 
Jo 

We define the geodesic distance d(x,y) as the infimum 

d(x, y) = mi Len7 (1-3) 
for all locally Lipschitz paths 7 connecting x to y. 

( 2 'lf M is infinite dimensional, then we suppose that it is modeled on a Banach or Hilbert 
separable space (see [Lan99], ch.II). 

( 3 ) Sometimes F is called a "Minkowsky norm" 

( 4 ) As suggested in [MM], we want to avoid referring to 7 as a curve, because confusion arises 
when we will introduce the manifold M of closed curves. So we will always talk of paths in 
the infinite dimensional manifold M. Note also that these paths are open-ended, while curves 
comprising M are closed. 
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The path connecting x to y that provides the min.E(7) in the class of all 
paths 7 connecting x to y is the (minimal length) geodesic connecting x to y. 

Note that, in the classical books on Finsler Geometry (see for example 
[BCS]), M is finite dimensional, and F(c, v) 2 is considered to be smooth and 
strongly convex in the v variable (for v ^ 0) (see also 1.1.2); this hypothesis is 
not needed, though, for the theorems in §A. 

1.1.1 Riemannian geometry 

Definition 1.4 Suppose that M is modeled on a Hilbert separable space H. A 
Riemannian geometry is defined by associating to M a Riemannian metric g; 
for any c G M , g(c) is a positive definite bilinear form on the tangent space 
T C M of M at c: that is, if h, k are tangent to M at c, then g(c) defines a 
scalar product (h,k) c . If the form g is positive semi-definite then the geometry 
is degenerate, and we will speak of a pseudo-Riemannian metric. If it is positive 
definite, then tangent space T X M is isomorphic to H , by means of the metric g. 
See [Lan99], ch. VII. 

A Riemannian geometry is a special case of a Finsler geometry: we define the 
norm 

\h\ c = ^thJT) , 

(pseudo-Riemannian geometries produce a seminorm \h\ c ). 

In the following we often drop the base point c from (h, k) c and \v\ c . 

If M is finite dimensional, then we can write 

(h,k) c = hig i>j (c)kj 

in a choice of local coordinates; then the matrix g 1,3 (c) is smooth and positive 
definite. 

1.1.2 Geodesies and the exponential map 

Suppose that the metric F is a regular Finsler metric, that is, F is of class C 2 
and F(c, -) 2 is strongly convex (for v ^ 0); such is the case when F(x,v) 2 = 
\v\ 2 = (v,v) in a smooth Riemannian manifold. 
Let 

i = r(7,7) 

be the Euler-Lagrangc O.D.E. characterizing critical paths 7 of 

E(i) - f \j(v)\ 2 dv 
Jo 

Define the exponential map exp c : T p M — > M as exp c (ry) = 7(1) when 7 is 
the geodesic curve solving 

{7(i>) =r(7(«), 7(«)), 7(0) = c, j(0)=V (1.5) 
Then we may state this extended version of the Hopf-Rinow theorem 

( 5 )As explained in the proposition A.l in the appendix A, we may equivalently define a 
minimal geodesic to be a minimum of minLen(7), that has been reparameterized to constant 
velocity 
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Theorem 1.6 (Hopf-Rinow) Suppose that M is finite dimensional and con- 
nected, then these are equivalent: 

• (M, d) is complete 

• closed bounded sets are compact 

• M is gcodcsically complete at a point c, that is, (1.5) can be solved for 
all v € IR and all r\, that is, the map rj i— > exp c (77) is well defined 

• for any c, the map rj i— » exp c (?y) is well defined and surjective; 

and all those imply that Vx, y G M £/iere exisi a minimal geodesic connecting x 
to y. 

1.1.3 Submanifolds 

The simplest examples of Ricmannian manifolds are the submanifolds of a 
Hilbcrt space H. We think of the finite dimensional case, when H = R n , 
or the infinite dimensional case, when wc assume that H is separable. 

Proposition 1.7 Define the distance d{x,y) = \\x — y\\ in H. Suppose that 
M C H is a closed submanifold. We may view M as a metric space, (M, d) : 
then it is complete. 

We may moreover induce a Riemannian structure on M using the scalar 
product of H: this in turn induces the geodesic distance d 9 , as defined in (IS). 
Then d 9 > d, and (M,d 9 ) is complete as well. If M is of class C 2 , moreover, 
then d and d 9 are locally equivalent. ^ 

It is not guaranteed that d and d 9 are globally equivalent, as shown by this 
example 

Example 1.8 (7 *> Let H = R 2 and M = {(s, sin(s 2 ))}. Let x n = {y/¥a,Q) <E 
M. Then d(x n , x n +i) — ► whereas d 9 (x n , x n +i) > 2. 

In a certain sense, infinite dimensional Ricmannian manifolds are simpler 
than their corresponding finite-dimensional: indeed, by [EE70], 

Theorem 1.9 (Eells-Elworthy) Any smooth differ entiable manifold M mod- 
eled on an infinite dimensional Hilbert space H may be embedded as an open 
subset of a Hilbert space. 

With respect to geodesies, the matter is though much more complicated. 
Suppose M is infinite dimensional. In this case, if (M, d) is complete the equa- 
tion (1.5) of geodesic curves can be solved for all v € IR; but (unfortunately) 
many other important implications contained in the Hopf-Rinow theorem are 
false. 

The most important example is due to Atkin [Atk75]: 

Example 1.10 (Atkin) There exists an infinite dimensional complete con- 
nected smooth Riemannian manifold M and x,y € M such that there is no 
geodesic connecting x to y. 

Proof by standard arguments, see for example sec. VIII. §6 in [Lan99] 
( 7 )\Ve thank A.Abbondandolo for this remark. 
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We show a simpler example, of an infinite dimensional Riemannian manifold 
M such that the metric space (M, d) is complete, but there exist two points 
x, y <E M that cannot be connected by a minimal geodesic. 

Example 1.11 (Grossman) Letl 2 (¥$) be the Hilbert space of real sequences 
x = (xq, X\, . . .) with the scalar product 



Let N = e = (1,0,0,.. .), S = -e = (-1,0,0,...). 

Let 7j be the geodesic starting from 7i(0) = N and with starting speed 7j(0) = 
e,; then there exists a first Xi > swc/i i/ia£ 7i(A,) = 5 (moreover Len(ji) = A J. 

TTien Len(^i) — > 7r 7 owi t/ie sequence 7, rfoes noi /lave a /imii. 

TVote i/iai me may iteifc 0/ using weak convergence: but 7^ weakly converges 
to the diameter; and weakly converges to 0. 

See also [Eke78]. 

It is then, in general, quite difficult to prove that an infinite dimensional 
manifold admits minimal geodesies (even when it is known to be complete); a 
known result is 

Theorem 1.12 (Cartan-Hadamard) Suppose that M is connected, simply 
connected and has seminegative curvature; then these are equivalent: 

• (M, d) is complete 

• for a c, the map i] —> exjp c (rj) is well defined 

and then there exists an unique minimal geodesic connecting any two points. ( 9 ' 

1.2 Geometries of curves 

Now, suppose that c(6) is an immersed curve c : Si — * R™, where S 1 is the 
circle; we want to define a geometry on M, the space of all such immersions c. 

M is a manifold; the tangent space T C M of M at c contains all the defor- 
mations h E T C M of the curve c, that are all the vector fields along c. Then, 
an infinitesimal deformation of the curve c in "direction" h will yield (on first 
order) the curve c(u) + eh{u). 

If 7 : [0, 1] — > M is a path connecting curves, then we may define a homotopy 
C : S 1 x [0, 1] -> ET associated to 7 by C(0, u) = j(v)(0) (more is in §3.2.1). 

( 8 'This example is also in a remark in sec. VIII. §6 in [Lan99] 
(^Corollary 3.9 and 3.11 in sec. IX.§3 in [Lan99] 



00 




Let e, = (0, 0, . . . , 0, 1,0,.. .) where 1 is in the i-th position. 
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1.2.1 Finsler geometry of curves 

Any energy that we will study in this paper can be reconducted to this general 
form 

E{ 1 ) = J ^( 7 (., w ),d„ 7 (., w )) 2 dv (1.13) 

where F(c, h) is defined when c is a curve, and h G T C M is a deformation of c; 
note that F will be often a Minkowsky seminorm and not a norm ' 10 ' on the 
space M of immersions (see 1.18). 

We look mainly for metrics in the space M that are independent on the 
parameterization of the curves c: to this end, [MM] define 

Bi = Bi(S\]R 2 ) = Imm(5' 1 ,R 2 )/Diff(5' 1 ) 

and 

Bij = BijiS 1 ,^ 2 ) = Imm / (5' 1 ,R 2 )/Diff(5 1 ) 

that are the quotients of the spaces Immj of smooth immersion, and of the 
space Imm/(S' 1 , IR 2 ) of smooth free immersion, with Diff(S' 1 ) (the space of au- 
tomorphisms of S 1 ). Bi j is a manifold, the base of a principal fiber bundle, as 
proved in §2.4.3 in [MM], while Bi is not. Any metric that does not depend on 
the parameterization of the curves c (as defined in eq.(1.15)) may be projected 
to Bi by means of the results in §2.5 in [MM] (the most important step appears 
also here as 3.10). 

Remark 1.14 (extending M) Imm(S 1 ,M) is on open subset of the Banach 
space C 1 (S' 1 — > IR"), where it is connected iff n > 3; whereas in the case n = 2 
of planar curves, it is divided in connected components each containing curves 
with the same winding number. 

To define a Riemannian Geometry on M — Imm{S 1 ,M), it may be conve- 
nient to view it as a subset of an Hilbert space such as iJ 1 (S' 1 — » R n ), and to 
complete it there. 

1.3 Abstract approach 

As a first part of this paper, we want to cast the problem in an abstract setting. 
There are some general properties that we may ask of a metric defined as in 
sec. 1.2.1. We start with a fundamental property (that is a prerequisite to most 
of the others). 

0. [well-posedness and existence of minimal geodesies] The Finsler 
metric F induces a well defined geodesic distance d; (M, d) is complete 
(or, it may be completed inside the space of mappings c : S 1 — > R n ); for 
any two curves in M, there exists a minimal geodesic connecting them. 

Let C be a minimal geodesic connecting cq and Ci. We assume that C is a 
homotopy of cq to c\ . 

( 10 JThat is, it will fail to satisfy property 3 in definition 1.1 

( 11 )Xhat is, the distance between different points is positive, and d generates the same topol- 
ogy that the atlas of the manifold M induces 
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1. [rescaling] if A > 0, if we rescale cq,ci to Xcq,Xci, then we would like 
that AC be a minimal geodesic. 

A sufficient condition is that F(Xc, h) = X a F(c, h) for some a > 0, and all 
A>0. < 12 > 

2. [euclidean invariance] If we apply an euclidean transformation A to cq, 
ci , we would like that AC be a minimal geodesic connecting Acq to Ac\ 

If F(Ac, Ah) = F(c, h) for all c, h, then the above is satisfied. 

3. [parameterization invariance] there are two version of this: we define 

curve-wise parameterization invariance when the metric does not 
depend on the parameterization of the curve, that is 

F(c,h) =F(c,h) (1.15) 

when c(t) = c(cp(t)) and h(t) = h(tp(t)) are reparameterizations of 
c, h 

homotopy-wise parameterization invariance Define 

C(0,v) = C(<p(0,v),v) 

where <p : S 1 x [0, 1] — > S 1 , <p G C 1 , ip(-, v) is a diffemorphism for all 
fixed v. We would like that, in this case, E(C) — E(C). 
If F can be written as 

F = F(c,TT N h) (1.16) 

(that is, F depends only on the normal part of the deformation) and 
if it satisfies (1.15), then, by proposition 3.8, E(C) = E(C). 

In both cases, the geometric structure we are building depends only on 
the embedding of the curves, and not on the parameterization. 

The above properties defined in 1,2 are valid for all examples that we will 
show; property 3 is satisfied for some of them. Property is possibly the most 
important argument. 

Definition 1.17 Any metric satisfying the above 1,2,3 is called a geometric 
metric. 

Note that 

Remark 1.18 If F satifies (1.16) then F(c,-) is necessarily a seminorm and 
not a norm^ 13 ^ on the space M : we should then talk of a pseudo-Ricmannian 
geometry of curves. The projection of F to the space Bi j may be nonetheless a 
norm. 

( 12 ) \Ve could ask F(Xc,h) = l(X)F(c, h) for some I : IR+ — > IR+ monotone increasing, with 
l(x) > for x > 0; but, if I is continuous, then l(x) = x a for some a > 

( 13 ) That is, it does not satisfy property 3 in the definition 1.1 
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These other properties would be very important in applications in Computer 
Vision. < 14 ) 

4. [finite projection] There should be a finite dimensional approximation 
of our metric, for purposes numerical computation. 

A sufficient condition is that the energy E(C) should be well defined and 
continuous with respect to a norm of a Sobolov space W k,p (I) (with k £ IN 
and p G [1, oo)): this would imply that we may approximate C by smooth 
functions C h and E{C h ) -> E{C). 

5. [embedding preserving] if cq and C\ are embedded, we would like C(-,v) 
to be embedded at all v 

6. [maximum principle] In the following, suppose that curves are embed- 
ded in R 2 , and write c C c' to mean that c' is contained in the bounded 
region of plane enclosed by c. ' 15 - ) 

If Co C c' , c\ C c[, then we would like that exists a minimal geodesic C 
connecting c' to c[ such that C'(-,v) C C(-,i>) for v € [0, 1]. 

This is but another version of the Maximum Principle; it would imply that 
the minimal geodesic is unique. It is an important prerequisite if we want 
to implement numerical algorithms by using Level Set Methods. 

7. [convexity preserving] if cq and ci are convex, we would like C(-,v) to 
be convex at all v 

8. [convex bounding] we would like that at all v, (the image of) the curve 
C(-,v) be contained in the convex envelope of the curves Cq,Ci 

9. [translation] If a € 1R™ ; if c\ — a + cq is a translation of Co, we would 
like that the uniform movement C{9, v) = Co(0) + va be a minimal geodesic 
from Co to c\ 

So we state the abstract problem: ( 16 ) 

Problem 1.19 Consider the space of curves M , and the family Q of all Rie- 
mannian (or, regular Finsler) Geometries F on M J 17 ^ 

Does there exist a metric F £ Q , satisfying the above properties 0,1,2,3? 

Consider metrics F £ Q that may be written in integral form 

F(c, h) = j f(c(s), d s c(s), . . . , dj,c(«), h(s), . . . , dl,h(s))ds 
what is the relationship between the degrees i,j and the properties in this section? 

l 14 'Some of these properties are much trickier: we do not know sufficient conditions that 
imply them 

( 15 ) By Jordan's closed curve theorem, any embedded closed curve in the plane divides the 
plane in two regions, one bounded and one unbounded 

( 16 ) Solving this problem in abstract would be comparable to what Shannon did for commu- 
nication theory, where in [Sha49] he asserted there would exist a code for communication on 
a noisy channel, without actually showing an efficient algorithm to compute it. 

is non empty: see [Lan99]. By using f.9, it would seem that there exists Riemannian 
metrics F on M such that (M, F) is geodesically complete; we did not carry on a detailed 
proof. 
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2 Examples, and different approaches and re- 
sults 

We now present some approaches and ideas that have been proposed to define 
a metric and a distance on the space of curves; we postpone exact definitions 
to section §3.2. 

2.1 Riemannian geometries of curves 

A Riemannian geometry is obtained by associating to T C M the scalar product 
of an Hilbert space H of squared integrable functions. We actually have many 
choices for the definition of H . 

2.1.1 Parametric (non-geometric) form of H° 

• We may define 

H = H^S 1 -> E") = L 2 {S 1 -> R") 
endowed with the scalar product 

(h,k) = [ (h(9),k(6))d6 (2.1) 

JSi 

for all h,k € T C M; this is a common choice in analysis and geometry 
texts ( 18 - ) ; however, the resulting metric is not invariant with respect to 
parameterization of curves (see (3) on page 8) and is therefore not geo- 
metric. 

Remark 2.2 This is the most common choice in numerical applications: 
each curve is numerically represented by a finite number m of sample 
points; thereby discretizing the geometry of curves to the geometry o/IR" m . 

Therefore (1.13) takes the form 

/ \\i{v)\\ 2 H odv=[ [ \d vl (6,v)\ 2 d6dv (2.3) 
Jo Jo Js 1 

that is defined when 7 G i? 1 ([0,l] — ► M). The energy of a homotopy is 
then 

E S (C) = J\d v C\ 2 

Definition 2.4 We define the space 

H lfi (l -> ]R n ) = H 1 ([0, 1] -* H^S 1 -* R")) 
and define the norm on the above space H 1 ' to be 

[ [ h(0,v)\ 2 + \d v7 (6,v)\ 2 d6dv (2.4.*) 
Jo Js 1 

( 18 )See ch 2.3 in [KH82], or the long list of references at the end of II§1 in [Lan99] 
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then 

Proposition 2.5 H > is the space of all finite energy homotopies, and 
the norm (2-4-*) above is equivalent to the energy (2.3) on families 7 with 
fixed end points (by prop. 3.14). 

E S (C) is strongly continuous in H ' and is convex, and E S (C) is coercive 
(proof by using 3.14). So E s has a very simple unique minimal geodesic, 
namely, the pointwise linear interpolation 

C*(9,v) = (l-v)c (d)+v Cl (9) 
• As a second choice, we may define in H the scalar product 

(h,k)= [ (h(s),k(s))ds= f (h{6),k(6))\c(6)\d6 (2.6) 
Jo JSi 

where I is the length of c, ds is the arc infinitesimal, and c = dgc. ( 19 ) 

This scalar product does not depend on the parameterization of the curve 
c; but the resulting metric is still not invariant with respect to rcparamc- 
terization of homotopies (see 3 in 8). 

By projecting this metric onto Bij and lifting it back to M (using 3.10), 
we then devise an appropriate geometric metric, as follows. 

2.1.2 Geometric form of H° 

We then propose the scalar product 

(h,k)= [ (n N h(s),n N k(s))ds (2.7) 
Jo 

where tin is the projection to the normal space N to the curve (see 3.2). From 
now on, when we speak of the H° metric, we will be implying this last definition. 

Note that we may equivalently define the scalar product as in (2.6), and only 
accept in T C M deformations that are orthogonal to the tangent of the curve. 
This would be akin to working in the quotient manifold Bij; or it may be 
viewed as a sub-Riemannian geometry on M itself. 

This geometric metric generates the energy 

E N (C)= Jjn N d v C\ 2 \C\dedv (2.8) 

which is invariant of reparameterizations of homotopies (see 3 in page 8). By 
proposition 3.10, the distance induced by this metric is equal to the distance 
induced by the previous one (2.6). 

( 19 'Thcre ; s an abuse of notation in (2.6); we would like, intuitively, to define h and k 
on the immersed curve c; to this end, we define h(s) and k(s) on the arc-parameterization 
c(s) : [0,Z] — > E" (with I = len(c)), and pull them back to c(0) : S 1 -» E", where we write 
h(0), k{0) instead of h(s(0)), k{s{0)) 
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Unfortunately, it has been noted in [Mum] that the metric H° does not 
define a distance between curves, since 

m£E N (C) = 

(see C.l and C.3). We will study this metric further in section 4. 

There is a good reason to focus our attention on the properties of this metric 
(2.7) for curves. Namely, this is precisely the metric that is implicitly assumed 
in formulating gradient flows of contour based energy functionals in the vast 
literature on shape optimization. Consider for example the well known geometric 
heat flow (dtC = d ss C) ( 20 - 1 in which a curve evolves along the inward normal 
with speed equal to its signed curvature. This flow is widely considered to be the 
gradient descent of the Euclidean arclength functional. Its smoothing properties 
have led to its widespread use within the fields of computer vision and image 
processing. The only sense, however, in which this is truly a gradient flow is 
with respect to the H° metric as we see in the following calculation (where L{t) 
denotes the time varying arclength of an evolving curve C(u,t) with parameter 

ue [0,1]). 

L(t)= I ds= ( \C u \du 
Jc Jo 

L'(t)= \ 1( ^-^du= Cc tu -C s du 

JO \Cu\ Jo 
\ / H° 

If we were to change the metric then the inner-product shown above would 
no longer correspond to the inner-product associated to the metric. As a conse- 
quence, the above flow could no longer be considered as the gradient flow with 
respect to the new metric. In other words, the gradient flow would be different. 
We will consider this consequence more at length in §5. 

2.1.3 Michor-Mumford 

To overcome the pathologies of the H° metric, Michor and Mumford [MM] 
propose the metric 

G*(h,k)= [ {l + A\ Kc \ 2 )(ir N h(u),ir N k(u))\c(u)\du (2.9) 
Jo 

on planar curves, where k c (u) is the curvature of c(u), and A > is fixed. This 
may be generalized to the energy' 21 -' 

E A (C)= f f \Tr N d v C\ 2 (l + A\H\ 2 )\C\ d6dv = E N {C) +AJ(C) (2.10) 
Jo Js 1 

on space homotopies C(u, v) : [0, 1] x [0, 1] — ► E™. This approach is discussed 
in §4.4.1. 

( 2 0)\ye will sometimes write C v for d v C (and so on), to simplify the derivations 
( 21 > J(C) is defined in (4.2); H is the mean curvature, defined in 3.3. Note that both E N (C) 
and J(C) are invariant with respect to reparamctcrization, in the sense defined in 3 on page 8. 



C t ■ C su du 



C t ■ C ss ds 



2.1 Ricmannian geometries of curves 
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2.1.4 Srivastava et al. 

We consider here planar curves of length 2ir and parameterized by arclength, 
using the notation £(s) : S 1 -> TEL 2 . Such curves are Lipschitz continuous. 
If £ £ C 1 , then |£| = 1, so we may lift the equality 

i(s) = (cos(0(«)), sin(0( S ))) (2.11) 

to obtain a continuous function 9 : 1R — > R. This continuous lifting is unique 
up to addition of a constant 2wh, h £ Z; and (9(s + 27r) — 9(s) = 2iri, where 
i £ Z is the winding number, or rotation index of £. The addition of a generic 
constant to 9 is equivalent to a rotation of £. We then understand that we 
may represent arc-parameterized curves £(s), up to translation, scaling, and 
rotations, by considering a suitable class of liftings 9(s) for s £ [0, 2-71"]. 

Two spaces are defined in [KSMJ03]; we present the case of "Shape Repre- 
sentation using Direction Functions'' ', where the space of (pre)-shapcs is defined 
as the closed subset M of L 2 = L 2 ([0, 2tt]), 

M = {9£ L 2 ([0,2tt]) I 0(0) = (2tt 2 ,0,0)} 
where : L? — > R 3 is defined by 

/•27T />27T /•27T 

01 (0) = / 9(s)ds , 02 (f) = / cos9(s)ds , 03 (6>) = / sin0(s)ds 
Jo Jo Jo 

Define Z as the set of representations in M of flat curves; then Z is closed 
in L 2 , and M \ Z is a manifold: ' 22 ^ 

Proposition 2.12 By the implicit function theorem, M \ Z is a smooth im- 
mersed submanifold of codimension 3 in L 2 . 

Note that M\Z contains the (representation of) all smooth immersed curves. 

The manifold M\Z inherits a Ricmannian structure, induced by the scalar 
product of L ; geodesies may be prolonged smoothly as long as they do not meet 
Z . Even if M may not be a manifold at Z, we may define the geodesic distance 
d 9 (x,y) in M as the infimum of the length of Lipschitz paths 7 : [0, 1] — ► L 2 
whose image is contained in M; ( 23 ^ since d 9 (x, y) > \\x — ?/||l 2 , and M is closed 
in L 2 , then the metric space (M, d 9 ) is complete. 

We don't know if (M, d 9 ) admits minimal geodesies, or if it falls in the 
category of examples such as 1.11. 

For any 9 £ M, it is possible to reconstruct the curve by integrating 

£(s) = [ cos(9(t)),sm(9(t)))dt (2.13) 
Jo 

This means that 9 identifies an unique curve (of length 2n, and arc-parametcrized) 
up to rotations and translations, and to the choice of the base point £(0); for 
this last reason, M is called in [KSMJ03] a preshape space. The shape space S 

( 22 ) The details and the proof of 2.12 and 2.15 are in appendix §B 

( 23 ) lt seems that M is Lipschitz-arc-connected, so d 3 (x,y) < 00; but we did not carry on a 
detailed proof 
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2 EXAMPLES, AND DIFFERENT APPROACHES AND RESULTS 



is obtained by quotienting M with the relation 9 ~ 9 iff 9(s) = 9(s — a) + b, 
a,6eE. We do not discuss this quotient here. 

Wc may represent any Lipschitz closed arc-parameterized curve £ using a 
measurable 9 € M: let arc : S 1 — * [0, 2n) be the inverse of 6 i— > (cos(0), sin(0)); 
arc() is a Borel function; then ((arc o £)(s) + a) G M, for an a e R. We 
remark that the measurable representation is never unique: for any measurable 
4,BC [0,27r] with |A| = \B\, (0(s) + 2n(l A (s) - l s (s))) will as well represent 
£ in M. This implies that the family A^ of that represent the same curve £ is 
infinite. It may be then advisable to define a quotient distance as follows: 

^=e^ A /^ ^ 

where d(9,9') = \\9 — 9'\\ L 2, or alternatively d = d 9 is the geodesic distance 
on M. 

If £ € C 1 , we have an unique^ 24 ' continuous representation 9 £ M; but note 
that, even if g C 1 , the infimum (2.14) may not be given by the continuous 
representations 9,9' of Moreover there are curves £ that do not admit a 

continuous representation 9. As a consequence, it will not be possible to define 
the rotation index of such curves £; indeed we prove this result: 

Proposition 2.15 For any h G Z, the set of closed smooth curves £ with rota- 
tion index h, when represented in M using (2.11), is dense in M \ Z. 

2.1.5 Higher order Riemannian geometry 

If we want an higher order model, we may define a metric mimicking the defi- 
nition of the Hubert space H 1 , by defining 

(h, k) H i = (h, k) H a + (h, k) H o (2.16) 

We have again many different choices, since 



we may use in the RHS of (2.16) the parametric H° scalar product (2.1), 
in which case the scalar product (h,k) H i is the standard scalar product 
of H 1 ^ 1 -> IT); then homotopics are in the space 

^([0, 1] — > H 1 (S 1 — > R")) 

with norm 



• we may use in the RHS of (2.16) the scalar product (2.6) or (2.7). Unfortu- 
nately none of these choices is invariant with respect to rcparamcterization 
of homotopics. 

We don't know of many application of this idea; the only one may be con- 
sidered to be [You98]. 



( 24 )lndeed, the continuous lifting is unique up to addition of a constant to 9(s), which is 
equivalent to a rotation of £; and the constant is decided by <j>i(9) = 2n 2 



2.2 Finslcr geometries of curves 
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2.2 Finsler geometries of curves 

To conclude, we present two examples of Finsler geometries of curves that have 
been used (sometimes covertly) in the literature. 

2.2.1 L°° and Hausdorff metric 

If we wish to define a norm on T C M that is modeled on the norm of the Banach 
space i°°(S' 1 -> ET), wc define 

F°°(c, h) = ||7rjv/i||z,°c, = supess e \-K^h{6)\ 

This Finsler metric is geometric. The length of a homotopy is then 

Len(C) = / supess e \-kn9 v C{0, v)\dv 
Jo 

Hausdorff metric We recall the definition of the Hausdorff metric 
du (A, B) = max < ma.xd(x, B), m&xd(A, y) > 

where A, B £ R n arc closed, and 

d(x, A)= min \x — y\ 

Let 3 be the collection of all compact subsets of 1R™. We define the length 
of any continuous path £ : [0, 1] — > S by using the total variation, as follows 

j 

lcn ff 7 = sup V d H (f (2.17) 
T 7=1 

where the sup is carried out over all finite subsets T = {to, ■ ■ ■ , tj} of [0, 1] and 
t < ■■ ■ < tj. 

The metric space (S, djj) is complete and path-metric, ^ 25 ^ and it is possible 
to connect any two A, B £ S by a minimal geodesic, of length dji{A, B). ( 26 ) 

Let S c be the class of compact connected A c R™; S c is a closed subset of 
(3, djj); 5c is Lipschitz-path-connected *- 27 ^; for all above reasons, it is possible 
to connect any two A, B £ 5 C by a minimal geodesic moving in 3 C ; but note that 
5 C is not geodesically convex in 5. ( 28 ) We don't know if (3 C , dn ) is path-metric. 



Projection of F°° into Hausdorff metric space When we associate to 
a continuous curve c £ M its image Jm(c) C R™, we are actually defining a 
natural projection 

Im : M — > 3 C 

f 25 ^ Path-metric: <1h(A,B) = inflen^7 where the infimum is computed in the class of Lip 
curves 7 : [0, 1] — * E connecting A to B 

(26) To prove this, note that is locally compact and complete; and apply A. 2 

( 27 ) That is, any A,B £ H c can be connected by a Lipschitz arc 7 : [0, 1] — > S c 

( 2S ) There exist two points A, B £ H c and a minimal geodesies J connecting A to B in the 
metric space (H, d), such that the image of £ is not contained inside H c 
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This projection transforms a path 7 in M into a path £ : [0, 1] — * S c ; if the 
homotopy C(9,v) = j(v)(6) is continuous, then £ is continuous. Moreover the 
projection of all embedded curves c is dense in S c . 
It is possible to prove that 

Len( 7 ) =len ff (£) 

for a large class of paths; and then the distance induced by the metric F°° 
coincides with df/(im(co), Im(c\)). 

This is quite useful, since in the metric space generated by the metric F°°, 
it is possible to find two curves that cannot be connected by a geodesic (this is 
due to topological obstructions); whereas, the minimal geodesic will exist in the 
space (S c , d H ). 

For this reason, [CFK03] proposed an approximation method to compute 
len ff (£) by means of a family of energies defined using a smooth integrand ^ 29 - ) , 
and successively to find approximation of geodesies. 

Unfortunately, the geometry in (S, dn) is highly non regular: for example, 
it is possible to find two compact sets such that there arc uncountably many 
minimal geodesies connecting them. ( 30 ' 



2.2.2 L 1 and Plateau problem 

If we wish to define a geometric norm on T C M that is modeled on the norm of 
the Banach space L 1 (S' 1 — ► R"), we may define the metric 

F\c,h) = \\n N h\\ L , = J \n N h(6)\\c(6)\d6 

The length of a homotopy is then 

Len(C) = J \ir N d v C(e, v)\\C(9, v)\d9dv 

which coincides with 

Len(C) = J \d v C(0,v) x d e C(0,v)\d8dv 

This last is easily recognizable as the surface area of the homotopy (up to 
multiplicity); the problem of finding a minimal geodesic connecting cq and c\ in 
the F 1 metric may be reconducted to the Plateau problem of finding a surface 
which is an immersion of / = S 1 x [0, 1] and which has fixed borders to the 
curves Cq and C\. The Plateau problem is a wide and well studied subject upon 
which Fomenko expounds in the monograph [Fom90]. 



3 Basics 

3.1 Relaxation of functionals 

To prove that a Riemannian manifold admits minimal geodesies, we will study 
the energy £^(7) by means of methods in Calculus of Variations; we review some 
basic ideas. 

( 29 'The approximation is mainly based on the property H/H^p —> p , for any measurable 

function / defined on a bounded domain 
( 30 'This is an unpublished result due to Alcssandro Duci. 



3.2 Curves and notations 
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Let X be a topological space, endowed with a topology r, and OcX, and 
/ : Q -> E. 

The function / is lower semi continuous if, for all x G X, 

liminf /(y) > /(x) 

J,— >x 

Note that, if the topology r is not defined by a metric, then we can introduce 
a different condition: F is sequentially lower semi continuous if, for all x <E X, 
for all x n ^ x, 

liminf /(x„) > /(x) 

n 

We define the sequential relaxation Tf of f on $1 to be the function 

17 : H 

that is the supremum of all /' : SI — * ]R that are sequentially lower semicontin- 
uous on SI, and f '< f in SI. We have that, for all x G 0, 

Tn/(x) = min {liminf /(x„)} 

(x n ).x n — >x n 

where the minimum is taken on all sequences x n converging to x, with x n , x G SI. 
f is sequentially lower semicontinuous, iff r/|n = /• 

If X is a metric space, then "sequentially lower semicontinuous functions" are 
"lower semicontinuous functions" , and viceversa; so we may drop the adjective 
"sequentially" . 

Consider again a function / : X — > R: it is called coercive in X if VM € IR, 

/(*) < M} 

is contained in a compact set. 

Proposition 3.1 If f is coercive in X and sequentally lower semi continuous 
then it admits a minimum on any closed set, that is, for all C C X closed there 
exists x G C s.t. 

f(x) = min f(y) 

This result is one of the pillars of the modern Calculus of Variations. We will 
see that unfortunately this result may not be applied directly to the problem at 
hand (see 4.1). 

3.2 Curves and notations 

Consider in the following a curve c(6) defined as c : S 1 — > H™. 
We write c for dgc. 

Definition 3.2 Suppose that c is C 1 ; or suppose that c G W,' , that is, c admits 
a weak derivative c. At all points where c{9) ^ 0, we define the tangent vector 

T(S\ t{9) 

At the points where c — we define T = 0. 
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Let v £ ]R™ . We define the projection onto the normal space N = T 

TT N V = v — (v, T)T 

and on the tangent 

tttv = (v, T)T 

so TrpfV + tttv = v and \itnv\ 2 + |tttv| 2 = \v\ 2 (that implies \ttnv1 2 = \v\ 2 — 

If c admits the weak derivative dec then T is measurable, so T G L°° and 
7rr,7r i vGi 00 (5 1 ^IR™ x ") 

A curve c £ C 1 (5 1 — ► R") is immersed when |c| > at all points. In this 
case, we can always define the arc parameter s so that 

ds = \c{9)\d6 

and the derivation with respect to the arc parameter 

We will also consider curves c £ W 11 (S 1 — » R n ) such that |c| > at almost 
all points. 

There are two different definitions of curvature of an immersed curve: mean 
curvature H and signed curvature k, which is defined when c is valued in M 2 . 

H and k are extrinsic curvatures ( - 31 ^: they are properties of the embedding 
of c into nr. 

Definition 3.3 (H) If c is C 2 regular and immersed, we can define the mean 
curvature H of c as 

H = d s T= ±-d e T 
|c| 

In general, we will say that a curve c £ 

w ioA Sl ) admits mean curvature in 
the measure sense if there exists a a vector valued measure H on S 1 such that 

T{s)d s cj)(s) ds = - [ <j){s)H(ds) V0 G C^'iS 1 ) 
J i J i 

that is 

r T{0)d e <i>{0) d9 = - cj)(9) \c\H(d8) V<£ £ C 00 ^ 1 ) 

Note that the two definitions are related, since when c £ C 2 , the measure is 
H = d s T ■ ds. See also [Sim8S], §7 and §16. 

We can then define the projection onto the curvature vector H by 

*hv = —^{v,H)H 



( 31 )See also: Eric W. Weisstein. "Extrinsic Curvature." From MathWorld-A Wolfram Web 
Resource, http : //mathworld. wolf ram. com/ExtrinsicCurvature .html 



3.2 Curves and notations 
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Definition 3.4 (N) When the curve c is in R 2 , and is immersed, we can de- 
fine^ 32 ^ a unit vector N such that N _L T and N is n/2 degree anticlockwise 
with respect to T. In this case for any vector V £ R 2 , ttnV = N(N, V), and, 

\V 2 \=TT N V 

Definition 3.5 (n) if c is in R 2 then we can define a signed scalar curvature 

K = (H, N) 

If H is a measure, then k is a real valued measure defined by 

n „ 

k(A) = 5^ / Ni{6)dHi{de) . 
i=i Ja 

Note that \H\ = \k\. 

When c G C 2 (S' 1 — > R 2 ) is immersed, 

d s T = kN and d s N = -kT . 

Remark 3.6 When the curve c is in R 2 , and is immersed then 

(H, v) = kttnv 

whereas for immersed curves c in R" 

\(H,v)\<\H\\n N v\ (3.6.*) 

and we do not expect to have equality in general when n > 3, since H is only a 
vector in the (n — 1)- dimensional subspace N = T . 

3.2.1 Homotopies 

Let I = S 1 x [0, 1]. We define a homotopy to be a continuous function 

C(6,v) : S 1 x [0,1] -> R" 

This homotopy is a path connecting cq = C(-, 0) and c\ = C(-, 1) in the space M 
of curves: indeed any path 7 in M is associated to a homotopy C by C(8, v) = 

7(f)(0)- 

We extend the above definitions to homotopies C, isolating any curve c in the 
homotopy by defining c(6) — C{6, v) for the corresponding fixed v; for example, 
when C admits the weak derivative dgC then 

TT T ,TT N € ^(S 1 X [0,1] ^R™ X ") 

Remark 3.7 We extend the measure H(-,v) on S , (that is curvature of any 
single curve C(-,v)) to a Borel measure H on I, by 

H(A) = f H(A v )dv (3.7.*) 
Jo 

(32)rpj lere j g a gijgjjt abuse of notation here, since in the definition TV = T 1 - in 3.2 we defined 
AT to be a "vector space" and not a "vector" 
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where A v is the section of A. H can be defined equivalently using the formula 

J T(9,v)d e <t>(9,v) = - J (f>{9,v)\C(0,v)\H(d9,dv) V</> G C™{I) 

Moreover we define the length 

len(C)(u) = / \C(9,v)\ d9 
Js 1 

so that len(C) : [0, 1] — > R + is a function of v. 
3.3 Preliminary results 

Definition 3.8 (Reparameterization) We define the reparameterization of 
C to C as 

C(9,v) = C(ip(9,v),v) (3.8.*) 

where tp : S 1 x [0, 1] — > S 1 is C 1 regular with dgtp ^ 0. Then (by direct compu- 
tation) 

d e C(6, v) = dgC(r, v)d e tp(9, v) (3.8. * *) 

so that T = Tsign{dgip); and 

TTfjd v C(9, v) = n N d v C(T, v) (3.8.o) 

whereas 

n f d v C(6, v) = w T d v C(r, v) + C(r, v)d v ip(6, v) (3.8. o o) 

where r=(p(6, v). 

We may choose to reparameterize using the arclength parameter as in the 
following proposition. 

Proposition 3.9 (Arc parameter) For any C 1 regular homotopy C such that 
all curves 9 i— * C are immersed, there exists a ip as in (3.8.*) above such that 
\dgC\ is constant in 9 for any v (that is, there exists I : [0, 1] — > R such that 
\d e C(0,v)\=l(v)) 

Proof. We just choose 

<p(0,v) = ^ d J\C(t,v)\dt 

□ 

On the other hand, we may reparameterize to eliminate nTd v C 

Proposition 3.10 For any C 2 regular homotopy C such that all curves 9 i— > C 
are immersed there exists a ip as in (3.8.*) above such that TTj^d v C = 0. 



3.3 Preliminary results 
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Proof. Both tttOvC and Cd v ip are parallel to T: so, if d v (p = —(d v C,T)/\C\, 
then TTfd v C = 
The O.D.E. 

T=(p(e,v) (3.io.*) 

can be solved for v 6 [0, 1], since 

|C(t,«)| 

is periodic in r and continuous, and then is bounded: that is, 

max M < oo 
S 1 x[o,i] 

Defining ij) = d$ip = cp we compute 

} = d (a t ,c(r,«),r(r,t;)) 

<fo Vl lV} dOdv^ d6 \C(t,v)\ 

d f(d v C(T,v),T( T ,v)) 



drV |C(r,u)| 



so -0 solves 



I { Ur ic(r,«)| (3.10.**) 

[V(#,0) = l, 6> eS 11 

and then > at all times v. □ 

Note that ip is not unique: wc may change in (3.10.*) and simply require that 
<£>(■, 0) be a diffcomorphism. The above result is stated in section 2.8 in [MM]; 
there C is called a horizontal path, since it is the canonical parallel lifting of a 
path in Bit to a path in M. 

For example, consider the unit circle translating with unit speed along the 
x-axis, giving rise to the following homotopy C(9,v): 

C{9,v) = (i> + cos6>, sin 6) 

While this is certainly the simplest way to parameterize the homotopy, it does 
not yield a motion purely in the normal direction at each point along each curve. 
However, the following reparameterization of the homotopy yields exactly the 
same family of translating circles (and therefore the same homotopy) but such 
that each point on each circle along the homotopy flows exclusively along the 
normal to the corresponding circle (i.e. in the radial direction). 



C(0,v) = (x,y) 

(l-e 2l, ) + (l + e 2l, )cos<9 



x{6,v) 
y(0,v) 



(1 + e 2v ) + (1 - e 2v )cos9 
2e v sin 9 
(1 + e 2v ) + (1 - e 2 ")cos0 
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In figure 1 we see a comparision between the the trajectories of various points 
(fixed values of 6) along the translating circle in the original homotopy C and 
its reparamctcrization C. 




C 



Figure 1: Reparameterization to iTj,d v C = for the translation of the unit circle 
along the x-axis. The dotted line shows the trajectory of a point on the curve. 

The above result is suprising, kind of "black magic" : it seems to suggest 
that while modeling the motion of curves we can neglect the tangential part of 
the motion •Kr r d v C. Unfortunately, however, this is not the case. We begin by 
providing a simple example. 

Example 3.11 The curve C(-,v) is translating as in figure 2, with unit speed: 

c(e,v)=c {e) + ve 1 




D A 

Figure 2: Stretching the parameterization 

After the reparameterization 3.10, the points in the tracts BC and AD are 
motionless, whereas the parameterization in AB is stretching to produce new 
points for the curve, and in CD it is absorbing points: then, dgC — > oo if the 
curvature in AB goes to infinity. 

We now consider the math in more detail. 

Proposition 3.12 Suppose that n = 2 for simplicity: define the curvature k as 
per definition 3.5 on page 19. 

Suppose in particular that \C\ = 1 at all points: then 

d v \C\ 2 = = 2(d v d e C,T) 

By deriving (3.8.*) we obtain 

d g C(e,v)=deC(T,v)i;(e,v) 

where ip = dgip solves (3.10.**), that we can rewrite (using \C\ — 1) as 

d v il> = -ip(d v C, N}k 



3.4 Homotopy classes 
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(where C,K,N,T are evaluated at ((p(8,v),v)). 

Note that \ (d v C, N)\ = \Tr^d v C\. This implies that the parameterization of 
C will be highly affected at points where both k and iriydyC are big. 

So the above teaches us that there are two different approches to the repa- 
rameterization: we may either use it to control C or tvt9 v C, but not both. 

3.4 Homotopy classes 

3.4.1 Class C 

Given an energy and two continuous curves cq and ci, we will try to find a 
homotopy that minimizes this energy, searching in the class C so defined 

Definition 3.13 (class C) Let I = S l x [0, 1]. Let C be the class of all homo- 
topies C : / — * TEC 1 , continuous on I and locally Lipschitz in S 1 x (0,1), such 
that Co = C(-, 0) and c\ = C(-, 1). 

Such minimum will be a minimal geodesic connecting cq and c\ in the space of 
curves. In this class C we can state 

Proposition 3.14 (Poincarfe inequality) there are two constants a', a" > 
such that VCeCn H 1 ' , 

||C||!ri.o(j) ^ «' + a " J \W (3.14.*) 
a', a" > depend on cq and c\. 

3.4.2 Class F of prescribed-parameter curves 

In some of the following sections we will change our point of view with re- 
spect to the above assumption 3.13 by fixing a measurable positive function 
I : [0, 1] — > R + and restricting our attention to a family of homotopies such that 

\d e c{e,v)\ = i{v). 

Definition 3.15 (class F) The class F contains all of C € C such that 

\dgC(6,v)\=l(v) 

for all 9, v. 

In particular, len(C) = 2irl(v). 

Note that, by 4.3, if I is not Holder continuous then the class F will be 
too poor to be useful (for example, it will not contain smooth functions). Un- 
fortunately the class F is not closed with respect to the weak convergence in 

Proposition 3.16 Suppose that I is bounded and \dgCh(0,v)\ = l(v) for all 
h,6,v; assume that 1 < p < oo and dgCh — v V weakly in L P (I — > R n ), or 
p = oo and d e C n ^* V weakly-* in L°°(I -> Ft"). 
Then \V(0,v)\ < l(v). 

Proof. Since |<9gC/j| < sup^, we may always assume that dgCh — i * V weakly-* 
in L°°(I — > R ra ); the result follows immediatly from theorem 1.1 in [Dac82] □ 
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We conclude that 

Theorem 3.17 Let 1 < p < oo. The closure of F n W 1,p with respect to 
weak convergence W x,p is contained in the class F of all C G W 1,p such that 
C(-, 0) = Co and C(-, 1) = c\ are given and \dgC\ < 

Proof. Suppose (Ch) C F and — 1 C: we prove that CsF. (Ch) is bounded 
in W 1,p : by Rellich-Kondrachov theorem (see thm. IV. 16 [Bre86]) (Ch) is pre- 
compact in L p : up to a subsequence, we may assume that Ch — » C in L p and 
that C — > C on almost all points: then C(-, 0) = c and C(-, 1) = ci- □ 

Remark 3.18 FFe could just as well have defined a class of homotopies with 
prescribed lengths, where we fix a function I and include in the class all C such 
that len(C) = l(v) for all v. However, if the energy E is geometric, then using 
the reparameterization 3.9, we can always replace any such C with a C G F, 
and E(C) = E(C). 

3.4.3 Factoring out reparameterizations 

If we only consider curves c such that |c| = 1, then (as pointed in [MM]) the 
group of reparameterizations is S 1 k Z2, where 

• the group S 1 is associated to the change of the initial point c(0) for the 
parameterization, that is, the reparameterization r G S 1 changes the curve 
c to c(6 + t), 

• and Z2 means the operation of changing direction, that is, c(9) becomes 
c(-6) 

We extend the above to homotopies satisfying \d$C(9, v)\ = l(v); the second 
reparameterization is not significant; the first one does not affect the normal 
velocity tt n8 v C; so the reparameterization in the class F does not affect energies 
that depend only on ir^dvC (as per eq. (1.16)). 

4 Analysis of E N (C) 

In this section we will focus our attention on the energy 
E N (C)= J\TT N d v C\ 2 \C\d9dv 

which is associated with the geometric version of the H° metric. We will derive 
a result of existence of minimima of E N in a class of homotopies such that the 
curvature is bounded. 

Remark 4.1 (winding) Unfortunately, any energy that is independent of pa- 
rameterization cannot be coercive. It is then difficult to prove existence of 
geodesies by means of the standard procedure in the Calculus of Variations 3.1. 
Indeed, suppose that C is a homotopy, and define 



C k (0,v) = C(6 + k2irv,v) 



4.1 Knowledge base 
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Vfc G Z; then 

E N {C k )=E N {C) 

(As a special case, consider two curves cq = c±, and Ck(9,v) — cq(9 + k2irv) ; 
each Cfc is a minimal geodesic connecting Cq,Cx, since E N (Ck) = 0). 

This proves that E N is not coercive in H 1 ^), since J \d v Ck\ 2 — > oo when 
k — > oo. 



4.1 Knowledge base 

Some of the follwing results apply in the class F of homotopies with prescribed 
arc parameter \dgC(0, v)\ = l(v) while others apply in the more general class C. 

Consider the integral 

J(C)= J^\H\ 2 \n N d v C\ 2 \C\d8dv (4.2) 

where H is the curvature of C along 0. We now deduce the following result from 
[MM] 

Proposition 4.3 Suppose that the homotopy C is smooth and immersed: then 
by extending to R" the computation in sec. 3.3 of [MM] we get 



1 



— Vlen(C)(v) <- x j I (H,d v C)*\C\d0 



s 1 



and then, for any < v' < v" < 1, 



v/lcn(C)K)- y/lcn(C){v') < i jf (j (H, d v C) 2 \C\de\ ' dv < 



2 



v' JS 1 



1/2 



< V^ 77 ^ 7 / / (H,d v C) 2 \C\d8dv) < 



< V 7 ^) rj, 1 

< — v v — V 

2 

(by Cauchy-Schwarz and (3.6.*)) and this implies that i/len(C) is Holder con- 
tinuous when J(C) is finite. 

Proposition 4.4 (l.s.c. and polyconvex function) Let W, V G Ft™. Let 

W x V be the vector of all n(n — l)/2 determinants of all 2 by 2 minors of 
the matrix having W, V as columns. 

Consider a continuous function f : R™ x — > 1R such that 

f(W, V) = g(W, V, W x V) 

where g : ]R"( rl + 3 )' /2 — > ]f| j s convex: then f is polyconvex ( 33 '. 



( 33 )The more general definition and the properties may be found sec. 4.1 in [But89], or in 2.5 
and 5.4 in [Dac82] 
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Let p > 2, suppose f > 0: by theorem J^.1.5 and remark 4-1-6 in [But89] 



then 



f{d 9 C, d v C)d9dv 



is W 1 ^ -weakly-lower-semi-continuous. ( 34 ) This means that, if Ch — v C weakly 
in W 1 *, that is 

Ch^C, d v C h ^d v C dgCh^dgC (4.4.*) 

in LP, < 35 ) then 

liminf ^ /(^C hl 9„C\) > jf ^C) 

4.2 Compactness 

We now list some simple lemmas. 

Lemma 4.5 let C(9,v) — C{6 + tp(v),v) be a reparameterization; then 



inf / \7r T d v C\ 2 = 



Js 1 



Tr T d v C(e,v) - f TT T d v C(s,v) ds ) de dv (4.5.*) 

s 1 



Lemma 4.6 (Poincare inequality in S 1 ) If f : S 1 — > Ft™ is C 1 , i/ien 



max|/|-min|/|<± f \df\ 



so that 



For any a € R, a ^ 0, 



7'-/'*=/. 



\df\ 



max |/| -min|/| < - / |d/ + a| 



so that 



(4.6.*) 



(4.6.**) 



(4.6.o) 



(4.6. oo) 



Lemma 4.7 Suppose C <G C 2 . //we derive d v (\C\ 2 ) we get 

d v {\C\ 2 ) = 2(T,d v d B C)\C\ 

so,< 36 ) u>/ien |C*| ^ 0, 

^(|C|) = {T,d v d e C) 
If the curves are in C 2 and we derive dg(iTTd v C), we get 



(4.7.*) 



d g (ir T d v C) = d g (T,d v C) = (H,d v C)\C\ + (T,d e d v C) = (H,d v C)\C\ + d v (\C\) 

(4.7.**) 

(where (H,d v C) = kttn8 v C for curves in the plane). 

( 34 ) Sequentially weakly-* if p = oo 

( 35 ) We can write equivalently Cj, — ► C or C7, C in (4.4.*), thanks to Rellich-Kondrachov 
theorem (see thm. IV. 16 [Brc86]) 

(36) when \c\ = 0j then T = o hy 

our definition, so the equation (4.7.*) holds as well in the 

distributional sense 



4.3 Lower semicontinuity 
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Proposition 4.8 Let M > be a constant. Suppose that a C 2 homotopy 
C : I —> R" satisfies \C(6,v)\ = l(v) and 

(J (H,d v C)\C\ dO^j dv<M (4.8.*) 

where H is the curvature of C along 9. Then there exists a suitable reparame- 
terizations 

C{6,v) =C(9 + <p(v),v) (4.8.**) 



such that 



J \n T d v C\ 2 < 2M 



Proof. Suppose that I € C 1 (the general proof being obtained by an approxi- 
mation argument). We summarize derivations in 4.7, 

d g {7r T d v c) = (H,d v c)\c\ + dj 

Applying eqns (4.5.*) and (4.6.oo) (with a = —d v l) we obtain 



inf / \ir T d v C\ 2 = f f ir T d v C(8,v) - -f ir T d v C(8,v) 
v J Jo Js 1 J s 1 



d6 dv < 



, dO 

s 1 

([ \d e (ir T d v C) + a\ d$] <[ (f \(H, d v C)\C\ + dj - dj\ do) dv 



o KJs 1 I Jo xJs 1 



□ 



Note that the reparamctrization (4.8.**) may be viewed as an unwinding when 
confronted with 4.1. 

Remark 4.9 1 Note that, if\C'(9,v)\ = l(v) then 

(H,d v C)\C\ do) dv< f (H,d v C) 2 l(v) 2 < (sup Z)J(C) 



10 VS 1 / J I 

So a bound on J(C) should provide compactness. 

4.3 Lower semicontinuity 

Let a > 0, > 0, V, W G R", define 

e{W,V) = \n w ±V\ a \Wf 

and define 

E a ,p{C)= I e(d e C,d v C) 



Note that E N (C) is obtained by choosing a = 2, ft = 1. In general if ft = 1 then 
Ea,p{C) is a geometric energy (see 1.17). 

Let W x V be the vector of all n(n — l)/2 determinants of all 2 by 2 minors 
of the matrix having W, V as columns. The identity 

\tt w ^V\ 2 \W\ 2 = \V\ 2 \W\ 2 - (V, W) 2 = \Vx W\ 2 
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is easily checked ( 37 - ) . Let 

f{W,V) = \VxW\ . 
Note that / is a polyconvex function (see 4.4) and that 

e(W, V) = |7r vvi ^| Q |M/| /3 = |W x V| a |W| p_Q (4.10) 

and 

E a .MC) = J^(deC,d v C) a \d e Cf- a 

We can provide this lower semicontinuity result in the class F. 

Proposition 4.11 Let p > 2. Suppose a > j3 > 0. Fix a continous function 
I : [0, 1] — - » R + ; ^ > 0, and use it to build the class F. Then E a> p is W 1,p - 
weakly-lower-semi-continuous in the class F. 
More precisely: let 

C g F, (C*),, c F 
if Ch C weakly in W 1,p , that is, 

Ch — > C, 9„Ch — 1 (%C ^C/, — ^ <9gC 

in .L p ; and 

I(«) = |C(fl,i;)| = |C' fc (e,t;)| 

i/ien 

liminf E a>fl (C h ) > E a<f3 (C) 

h 

Proof. We prove the theorem in steps: 

• Let A > 0. Suppose \C\ = \C h \ = A; then 

E a ,fj(C) = \ f3 ~ a f f(d e C,d v C) a 



is l.s.c. (by 4.4). 

for any homotopy C and any continuous g : [0, 1] — > R + , define e g (C)(v) : 
[0,1]->]R+ 

„ I 9(vf- a / S i a„c)«^ if > o 

69l ° j "\0 if$(«) = 

If C e F (that is l(v) = \C(6,v)\) then 

E a ,p(C) = [ e^dv 
Jo 



( 37 )ln R 3 we have (V,W) = |V||W|cosa: and \V X W\ = |V||W|sina:, where a is the angle 
between the two vectors 
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Consider a piecewise function g > defined 

m 

9 = J2^X[a z ,a l + 1 ) (4.12) 



and let 

E(C)= e g (C)dv = y2 e gi (C)dv 

Jo i=1 Ai 

then we apply the previous reasoning to all addends and conclude that 

liminf E(C h ) > E{C) 

h 

• Suppose Z is continuous and I > 0. Let t be the class of piecewise functions 
g > defined as in (4.12), such that on any interval [a*, ctj+x), cither ^ = 
or 

.9; > sup Z . 

[(H,0«+l) 

Then for any such g and C G F, 



e g (C)dv < E aiP (C) . 

• choose g in the class r; then 

liminf E at p(Ch) > liminf / e g (Ch)dv > I e g {C)dv 



• Fix C: it is possible to find a sequence (pj) C t such that 

e w (C)(«) -^jei(C)(v) 
for almost all points v, monotonically increasing: indeed, let 

Aj,i = [i2- j ,(i + l)2- j ), I jti = inf l(v), Sj,i= sup l(v) 



and 

9j(v) = 

Then 



if v € Aj i and Ij^ = 
Sj i if v € A^j and Ij^ > 



E a ,0(C) = sup / e gj (C)dv 

3 JO 



We would like to prove this more general statement 



□ 
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Conjecture 4.13 choose Lipschitz homotopies 

c-.i^n n , c h -.i->n n 

define 

l(v)=lenC, lh(v)=lenCh 
If lh —> I uniformly and Ch — 1 C weakly in W l p , then 

limini E N (C h ) > E N (C) 

h 

Whereas we cannot generalize the theorem further: this is due to example 
4.3.1. 

4.3.1 Example 

We want to show that E N is not l.s.c. if we do not control the length. 
Let a > (3 > 0, define 

e(W,V) = \Ti w ±V\ a \Wf = \W x V\ a \Wf- a 

and define 

E a A G )= J<d e C,d v C) 
We will actually show that E a ^p is not l.s.c, and that 
Proposition 4.14 

TE a .p{C) = 

where the relaxation is computed with respect to weak- * L°° convergence of the 
derivatives. 

1. For simplicity, we temporarily drop the requirement that the curves C(-,v) 
be closed. We then redefine / = [0, 1] x [0, 1]. 

2. Let C : I — > / be a Lipshitz map such that 

(7(0,0) = (M), C(6, 1) = (6, 1), C(0,v) = (0,v), C(l, v) = (1, v) . 

(4.15) 

3. Let h > 1 be integers. 

We rescale C and glue many copies of it to build C\, as follows 

C ft (0,w)=ic , ((ft0)mod(l),(/i«)mod(l)) +b h (6,v) 
where bh '■ I — > H 2 is the piecewise continuous function 
b h (0,v)=(^[h0\,^[hv] 

(in particular, C\ = C). We represent this process in figure 3 on the facing 
page. 
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C c h 

Figure 3: Tesselation of homotopy C to form Ch 



Then 



d e C h (0, v) = d e C[{h6) mod(l), (hv) mod(l) 
d v C h (9, v) = d v c((hO) mod(l), (hv) mod(l) 



We may think of Ch as homotopies that connect the same two curves, 
namely, 

C h (6,0) = (6,0) = (7(0,0), C r fc (fl,l) = (fl,l)=e(fl,l) , 
while extrema points move in a controlled way, namely 

C h (0,«) = (0 J i;) = C(0 > t;) > C h (l, v) = (1, v) = (7(1, v) . 

4. Let C(6,v)=(9,v) be the identity. 

The sequence Ch has the following properties 

(a) Ch — > C in and more precisely 

sup\C h -C\ < \ 
i h 

(b) dgCh and 9t,C^ are bounded in L°°, and more precisely, 

sup \deCh\ < \deC\, sup \d v C h \ < sup 
/ / / 



and then all Ch are equi Lipschitz; 



(c) 



d e C h ^e 1 = (1,0) = d e C 



and 



d v C h ^e 2 = (0,l) = d v C 
weakly* in L°°(I). ^ 



( 3S ) Proof by lemma 1.2 in [Dac82] 
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(d) let l(v)=lea(C), a=J I and 

l h (v)=len(C h )(v) = l((hv) mod(l)) 

then l h a weakly* in L°°([0, 1]). < 38 ) 

(e) Suppose C is piecewise smooth, so that the curvature Hh of Ch can 
be defined almost everywhere. By (4.15), 1(0) = 1(1) = 1. If l(v) > 1 
at some points, then the sequence lh(v) is not equicontinuous: then, 
by proposition 4.3, the sequence of integrals 



(H h ,d v C h ) 2 \C h \d9dv 

o Js 1 



is unbounded in h. ^ 
(f) E a ^(Ch) is constant in h: indeed 

f i/h r i/h 

e(d e C h ,d v C h ) = h 2 / e(d e C h ,d v C h )d8dv = 

Jo Jo 
r i/h f i/h 

= h 2 / / e(d e C(h9,hv),d v C(h9,hv))d9dv 

IQ Jo 

i(d e C(9,v),d v C(9,v))d9dv 



that is 

E a ,p(C h ) = E a . p (C) (4.16) 
5. Let j > 1 be a fixed integer. Let C : / — > IR 2 be defined by 
7(«) = | " if ^ G [0, 1/2] 

71 j if »6 [1/2,1] { ^ U) 

and 

(5(0, u) = (0, t> + 7(v) sm(27rj(9)) (4.18) 

(note that C is a graph). 

Wc have this further properties 

• We know that lh — 1 a with a ~ J len C; but this limit a is strictly 
bigger than 1, whereas lenC = 1 

• For 9 G [0, l/h] and v G [0, l/h] 

d v C h (9,v)= d v C(h9,hv) = (0, 1 + j'(hv) sm(2-jr jh9)) 
d 9 C h (9,v)= d e C(h9,hv) = (1, 2irn(hv) cos(2njh9)) 



Then 



sup l^oCfal < 1 + 2irj, sup |^C h | < 2 



( 39 )a fortiori, J(C^) is unbounded — J{C) was defined in (4.2) 



4.4 Existence of minimal geodesies 
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• We compute 

E a ,p{&)= I e(d 9 C(9,v),d v C(9,v))d9dv 



o Jo 



|1 + 7» sin(j6»)| Q |l + j{v) 2 ] 2 coa(j9) 2 \^~ a ^ 2 d9 dv 



fl/2 pi 

2 I / \l + sm(j0)\ a \l + v 2 j 2 cos(j6) 2 \^- a V 2 d6dv 



then 

lim E af3 {C) = (4.19) 

J— *oo 

6. Combining (4.16) and (4.19) we prove that -E a ,/3 is not l.s.c. Indeed for j 
large 

Iim£? a|i8 (Cfc) - S Qj(3 (C) < E a .p{C) = 1 
whereas Ch — 1 C. 

7. to prove 4.14, consider any homotopy C; this may be approximated by a 
picccwise linear homotopy which in turn may be approximated by many 
replicas of the above construction. 

4.4 Existence of minimal geodesies 

Theorem 4.20 Let M > 0. 

Let A be the set of admissible curves c : S l — > R", suc/i i/iai 

• c : S* 1 — > R™ is Lipschitz, and c admits curvature H in the measure sense 
(see in 3.3), and moreover 

• the total mass l-f^KS* 1 ) of the curvature H of c is bounded uniformly 

\H\(S X )<M . (4.20.*) 

Let Co, ci be curves in A. 

Fix a bounded continuous function I : [0, 1] — ► R + , wi£/i inf Z > 0. 
Let B be the class of homotopies C : I — > R™ swc/i £/ia£ 

• C EH 1 ^ ^ R" ) 

• an?/ given curve 9 ^ C{9, v) is in A 

• i/ie curvature can be extended to the homothopy (see 3.7), d v C is contin- 
uous, and 



H,d v C)\C\d9) dv<M ((4.8.*)) 
Si / 

• l(v) = lenC(v) = J sl \C\ dO 

• C{9, 0) = c (9) and C(9, 1) = ci 

J/ $ is non empty, then the functional E admits a minimizing homotopy 
C* ; this minimum C* satisfies all the above reguirements, but possibly for con- 
dition (4-8.*). 
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Proof. The proof is divided in two important (and independent) steps 

• Let Ch be a sequence such that 

lhnE N (C h )=M B E N (C) 

Up to reparametcrization, assume 

\C h (8,v)\ = l(v) . 

By this bound, and the compactness result 4.8, we reparameterize any 
term Ch to Ch by 

C h (6,v) = C h (0 + (p h (v),v) 

so that 

J\7T T d v C h \ 2 <2M ; 

moreover 

J\7T N d v C h \ 2 < (maxy) J \n N d v C'h\ 2 \C\ < + M E N ) (4.21) 

(definitively in ft,) and then 

J \d v Ch\ 2 < 2M+(maxi)(l + inf£; Ar ) 

whereas ^ 

^|C7 ft | 2 = jf ;( V ) 2 d V : 

then (by the Banach-Alaoglu-Bourbaki theorem < - 40 - ) ) up to a subsequence, 
Ch converges weakly in H 1 to a homotopy C* . 

• We want to prove that \C*(8,v)\ = l(v) for almost all 8, v. We know that 
\C*(8,v)\ < l(v), by 3.16. Suppose on the opposite that \C*(8,v)\ < l(v): 
then there exists e > and an measurable subset A <Z I with positive 
measure, such that \C*(8,v)\ < l(v) — e for (8,v) e A. Let A v = {8 : 
(8, v) £ A} be the slice of A: then, by Fubini-Tonelli, there is a v such 
that the measure of A v is positive. We fix that v. Suppose that Hh is the 
curvature of Ch(-, v): then 

H h = l{v)dggC h 

in the sense of measures. We know that H has bounded mass: so dggCh 
has: by Theorem 3.23 in [AFPOO], dgCh(-,v) is compact in L 1 (S' 1 ), so, 
(up to a subsequence), we would have that Ch(-,v) — * C*(-) strongly in 
L 1 (5' 1 ): then \C*(8,v)\ = l(v) for that particular choice of v, achieving 
contradiction. 

• By 4.11, liminfh E N (C kh ) > E N (C*): then C* is the minimum. 

□ 

( 4 °)See thm 111.15 111.25 and cor 111.26 in [Bre86] 
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Remark 4.22 If we wish to extend the above theorem, we face some obstacles. 

• If we do not enforce some bounds on curvature (as (4-20.*) and (4-8.*)), 
then the example in %4-6 shows that we cannot achieve compactness of a 
minimizing subsequence 

• If we wish to remove the hypothesis "mi I > ", ' 41 ^ we are faced with 
the following problem: if l(v) = 0, then the curve C(-,v) collapses to a 
point; consequently if l(v) = on an interval [a,b], then the homotopy 
collapses to path on that interval; moreover, the length and the energy 
of C restricted to v £ [a, b] is necessarily 0, so that E(C) provides no 
bound on the behaviour of C : we again lose compactness of a minimizing 
subsequence (indeed, the inequality (4-21) needs that inf„ len(C)(u) > 0, 
to be able to control J \d v C\ ). 



4.4.1 Michor-Mumford 

Let d(co , c\ ) be the geodesic distance induced by the metric G A defined in 
[MM] (see eq.(2.9) here). By the results in [MM], we know that (in general) this 
distance is non degenerate: 

Proposition 4.23 Consider an homotopy C connecting two curves Cq = C(-, 0), C\ 
C(; 1), and its energy E A {C)=E N {C) + AJ(C) (see eq.(2.10) here) 
By Holder inequality, 

\C\] 1 (I \n N d v C\ 2 \C\) ' > [ |7r^C||C|> / \d v CxC\ 



We then obtain the area swept out bound ' 42 ^ 

\d v CxC\) 2 < EN{G) 



i 



Indeed the leftmost term is the area swept by the homotopy. 

By the proposition 4-3, we know that, if J(C) is bounded, then len(C) is 
continuous and is bounded. So, if Cq ^ C\, and there is no zero-area homotopy 
( 43 ) connecting Cq,C\, then d(ci,c 2 ) > 0. See also prop. 5.5 here. 

We cannot instead currently prove a theorem of existence of minimal geodesies 
for the energy E A (C); in this section we have though derived some insight; so 
we discuss the conjecture 

Conjecture 4.24 Fix two curves cq and c\. The energy E A (C) admits a min- 
imum in the class C of homotopies connecting cq and c\. 



May we improve the proof of theorem 4.20 to prove this conjecture? We 
discuss what is ok and what is wrong. 

( 41 ) note that the bound (4.20.*) does not imply that inf I > 

( 42 ) in a form slightly better than the one in §3.4 in [MM] 

( 43 ) note that there are curves co ^ ci in the completion of the space Bi (the completion is 
described in §3.11 in [MM]) that can be connected by a zero-area homotopy 



36 



4 ANALYSIS OF E N (C) 



• Wc may want to use J(C) to drop the requirement (4.20.*) that is, in 
turn, used to have l.s.c. of the functional E N ; so we may think of proving 
this lemma 

"Suppose that we are given a sequence of smooth homotopies Ch, and 
J(C h ) < M, and C h C weakly in W 1 *: then len(C h ) -► len(C)" 

but this is wrong, as seen by this example 

Example 4.25 Let C h : [0, 1] x [0, 1] -> R 2 defined as 

Ch(u,v) = ( u, — sin(27r/m) 

V h 

and 

C(u,v) = (u,0) . 

These homotopies do not depend on v: then J(Ch) = 0. On the other 
hand, Ch C but len(C^) is constant and bigger than 1 = len(C). 

• We need way to be sure that curves in the homotopy do not collapse to 
points (as discussed in 4.22); so we may think of proving this lemma 

"Suppose that the homotopy C admits curvature, and J(C) < oo: then 

inMen(C)(u) > 0" 

but this is wrong, as seen by this example 

Example 4.26 Let 

die) = (sin(0),cos(0)) 
be the circle in R 2 , and build the homotopy 

C(6,v) = v 4 Cl (9) 

Then 

J(C) 



-^(4v 3 ) 2 2irv i dv = 32tt^ v 2 dv = 32tt/3 



• We need a semicontinuity result on J(C). Indeed we cannot hope in any 
cancellation effect in the sum E N '(C) + J(C) because the two energies 
scale in different ways: ( - 44 ^ 

Remark 4.27 (rescaling) Lete>0 andC = eC thenE N {C) = e 3 E N (C) 
but J(C) = e J(C) 

• On the bright side, we do not have a counterexample to show if J(C) is 
not l.s.c; actually, remark (4e) (on page 32) suggests that if Ch — ► C and 
the homotopies have a common border condition (such as (4.15)), and 
J{C h ) is bounded, then liminf E N (C h ) > E N (C). 



( 44 'and this suggests that the energy E A (C) should not satisfy the rescaling property 1 on 
page 8 



4.5 Space of Curves 
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• Moreover, by 4.23 we know that the induced distance is in general not 
degenerate. 

• As pointed in 4.9, the term J(C) in the metric provides compactness in 
H\I) 

• Moreover the example in §4.6 shows that we do need to control the cur- 
vature of curves to be able to prove existence of minimal geodesies: this 
justifies the term J(C) in the Michor-Mumford energy E A (as well as the 
bound (4.8.*) in theorem 4.20). 

4.5 Space of Curves 

By using the previous theorem 4.20, we immediatly obtain a metric on a Space 
of Shapes. 

Fix M > 0. Let S be the space of unit length ( 45 ) closed immersed C 2 curves 
such that for any c e 5, the curvature k of c is bounded by M, as 

|k| < M . 

We may think of S as a "submanifold with border" in the manifold M of all 
closed unit length immersed C 2 curves. 
Then we can use the Ricmannian metric 

(h,k}= I (7r N h(e),Tr N k(9)) d6 

J Si 

to define a positive geodesic distance in S: by the theorem in the previous 
section, this distance admits minimal geodesies: 

Proof, indeed, we fix lcn(C) = 1; wc may write 

J(C) < K 2 E N (C) 

then we may use the remark 4.9 to obtain compactness of minimizing subse- 
quences. □ 

Unfortunately, since S has a border (given by the constraint |k| < M) then 
the minimal geodesic will not, in general, satisfy the Euler-Lagrange ODE de- 
fined by E N . 

4.6 The pulley 

We show that there exists a sequence of Lipschitz functions Ch : I — > R 2 such 
that |7TAf9„C^| < 1, \d e C h \ = 1, but fj \d v C h \ -> oo. 

Consider the figure 4 on the following page. The thick line ABCDEF is the 
curve i— ► C^{6, 0). The thick arrows represent the normal part ir^d-vCc, of the 
velocity, while the thin dashed arrows represent the tangent part -kt^vC^ of the 
velocity. The circles are just for fun, and represent the wheels of the pulley. 

( 45 )if the conjecture 4.13 holds true, then the "unit length" constraint may be dropped; note 
that the formula of the metric is "geometric" (as defined in 1.17); but the minimal geodesies 
and the distance are not invariant with respect to rescaling, due to the bound \k\ < M. 
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Figure 4: The pulley Ch, in case h = 5 



The movement of C§(9,v), that is, its evolution in v, is so described: the 
part ABCD is still, that is, it is constant in v; in the part DE (respectively 
FA) of the curve, vertical segments move apart (resp., together) as the thick 
arrows indicate, with horizontal velocity with norm lir^dyC^l = 1; as the curve 
unravels, it is forced to move also parallel to itself. 

The generic curve Ch has 2h wheels: h wheels in section DE, to pull apart, 
and h wheels in section AF, to pull together; the horizontal tracts in AF and 
DF are of lenght oc 1/h, so the tract AF straightens up in a time Av oc 1/h: 
at that moment, the movement inverts: so while v € [0, 1], the cycle repeats for 
h times. 

In this case, the tangent velocity in section DF is h times the normal velocity 
in sections AF and DE. Then, if we choose the normal velocity to have norm 
1, the tangent velocity will explode when h — ► oo. This means that the family 
of homotopies Ch will not be compact in H l {I — > R 1 ). 

The first objection that comes to mind when reading the above is "this ex- 
ample is not showing any problem with the curve itself, it is just giving problems 
with the parameterization of the curve" . 

Indeed we may reparameterize the curves so that the tangent velocity will 
not explode when h — > oo: by using 3.10, we obtain that wfd v C = 0. 

Remembering remark 3.12, we understand that this is not going to help, 
though. So we point out this other problem. 

Let A = X(v,h) be the distance from feature point D to feature point E. Then 

d v A — > oo if h —t oo. 

5 Conformal metrics 

We recall at this point that the well known geometric heat flow (C't = C ss ) < - 46 ' 1 
is truly the gradient descent for the Euclidean arclcngth of a curve with respect 
to the H° metric. Unfortunately, given the pathologies encountered thus far 
with H°, we see that this famous flow is not a gradient flow with respect to 
a well behaved Riemannian metric. If we propose a different metric, the new 

( 46 ) Recall that we write C v for d v C (and so on), to simplify the derivations 
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gradient descent flow for the Euclidean arclength of a curve will of course be 
entirely different. For example, the metric (2.9) proposed by [MM] yields the 
following gradient flow for arclength: 



Ct 



C s 



1 + A C ss ■ C s 



(5.1) 



Notice that the normal speed in (5.1) is not monotic in the curvature; and there- 
fore the flow (5.1) will not share the nice properties of the geometric heat flow 
(d t C = d ss C). For example, embedded curves do not always remain embedded 
under this new flow, as illustrated in Fig. 5. 





Figure 5: Intersections induced by flow (5.1): for big choices of A, the point x 
will travel faster than the point y (and in the same direction) and will eventually 
cross it. 

Given the pathologies of H° we have no choice but to propose a new metric 
if we wish to construct a well behaved Ricmannian geometry on the space of 
curves. However, we may seek a new metric whose gradient structure is as 
similar as possible to that of the H° metric. In particular, for any functional 
E : M — > M we may ask that the gradient flow of E with respect to our new 
metric be related to that gradient flow of E with respect to H° by only a time 
reparamctcrization. In other words, if C(t) represents a gradient flow trajectory 
according to H° and if C(t) represents the gradient flow trajectory according 
to our proposed new metric, then we wish that 



C(t) = C(f(t)) 

for some positive time reparameterization / : IR 
gradient flows will then be related as follows. 

Ct = f(t) C t 



R, / > 0. The resulting 



(5.2) 



The only class of new metrics that will satisfy (5.2) are conformal modifica- 
tions of the original H° metric, which we will denote by H®. Such metrics are 
completely defined by combining the original H° metric with a positive confor- 
mal factor <f> : M — > R where <f>(c) > may depend upon the curve c. The 
relationship between the inner products is given as follows. 



hi, h$ 



ii" 



0(c) (hi, h-i) 



(5.3) 
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Note that for any energy functional E of curves C(t) we have the following 
equivalent expressions, where the first and last expressions are by definition of 
the gradient and the middle expression comes from the definition (5.3) of a 
conformal metric. 

| £(c(t „ . (f ,v3£)\ - .v^c,) - (£. vbjq) 

\ Conformal ' H° \ Conformal 'if \ Original / JJ° 

Gradient Gradient Gradient 

(5.4) 

We see from (5.4) that the conformal gradient differs only in magnitude from 
the original H° gradient 

\/4> E = Ly E 

9 

and therefore the conformal gradient flow differs only in speed from the H° 
gradient flow. 

8C 1 

As such and as we desired, the solution differs only by a time reparameterization 
f given by 

7 0(C) 

The obvious question now is how to choose the conformal factor. 
A first suggestion is in the following theorem 

Proposition 5.5 ' 47 ^ Suppose that 

min(</>(c)/len(c)) = a > (5.5.*) 

c 

Consider an homotopy C connecting two curves cq = C(-,0),ci = C(-,l) ; and 
its H® -energy 

<f>(C(-,v)) I \ir N d v C\ 2 \C\d9dv 
Js 1 

Up to reparameterization 3.9, \C(d, v)\ = len(C(-, v))/2n so we can rewrite 
the energy (using (4-10)) as 

2 ?^[' ,V }} I \d v C xC\ 2 d0dv>2na f [ \d v C x C\ 2 d9 dv > 
o len(C(-,«)) Jgl Jo Jgl 

2 

~ a {j j \ d " C X C \ d0dyX 

the rightmost term is the square of the area swept by the homotopy. 

So, if Co 7^ c\, and there do not exist a homotopy connecting co,c\ with zero 
area, then d(c\,C2) > 0. 

Although we already know that the H° metric is not very useful, we may 
obtain a lot of insight into how to choose the conformal factor <p by observing 
the structure of the minimizing flow (which turns out to be unstable) for the 
H° energy in the space of homotopies. We will then try to choose the conformal 
factor in order to counteract the unstable elements of the H° flow. 



(47) "W e thanks prof. Mumford for suggesting this result. 



5.1 The Unstable H° Flow 
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5.1 The Unstable H° Flow 
5.1.1 Geometric parameters s and 

We have denoted by u € [0, 1] a parameter which traces out each curve in a 
parameterized homotopy C(u,v) and we have denoted by v € [0,1] the pa- 
rameter which moves us from curve to curve along the homotopy. Note that 
both of these parameters are arbitrary and not unique to the geometry of the 
curves comprising the homotopy. We now wish to construct more geometric 
parameters for the homotopy which will yield a more meaningful and intuitive 
expression for the minimizing flow we are about to derive. The most natural 
substitute for the curve parameter u is the arclength parmeter s. We must also 
address the parameter v, however. While v as a parameter ranging from to 1 
seems to have little to do with the arbritrary choice of the curve parameter u, 
the differential operator depends heavily upon this prior choice. The desired 
effect of differentiating along the homotopy is mixed with the undesired effect of 
differentiating along the contour if flowing along corresponding values of u be- 
tween curves in the homotopy requires some motion along the tangent direction. 
To see the dependence of on u, note that C(u, v) and C(u, v) where 



C(u,v) = c(u( 1+V \v^) 



constitute the same homotopy geometrically, and yet ^7 5^ frjf • 

We will therefore introduce the more geometric parameter whose corre- 
sponding differential operator yields the most efficient transport from one 
curve to another curve along the homotopy regardless of "correspondence" be- 
tween values of the curve parameters. It is clear that such a transport must 
always move in the normal direction to the underlying curve since tangential 
motion along any curve does not contribute to movement along the homotopy. 
More preceisely, we define the parameteres s and v* in terms of u and v as 
follows. 

d = 1 d md d = d ( C c) d 
ds \\C U \\ du dv* dv <9s 



5.1.2 H° Minimizing Flow 

Suppose we now consider a time varying family of homotopies C{u, v, t) : [0, 1] x 
[0, 1] x (0, oo) — » R" and compute the derivative of the H° energy along this 
family. Note that the H° energy, in terms of the new parameters s and u* may 
be simply expressed as follows (since ttnC v ^ = C Vt ). 

E{t) =11 \\C v ,fdsdv (5.6) 
Jo Jo 

In the appendix, we show that the derivative of E may be expressed as follows. 

E'{t) = -2 J J c t -(a^»-(a, w ,-c s )c s -(a ; c , ss )a,+^iia,ii 2 c ss )dsd ? ; 

(5.7) 

In the planar Cclsg, Cv. m cLnd Css are linearly dependent (as both are orthogonal 
to C s ) which means that 

{C Vt ■ C ss )C Vr = (Cv„ ■ C Vr )C ss = |jCu t |j C ss (5-8) 



42 



5 CONFORMAL METRICS 



and therefore 

E\t) = -2 J J C t - ((C v , v , - {C v , v , ■ C S )C S ) - \\\C v ,fC as ) dsdv (5.9) 
by which we derive the minimization flow 

Cf = C'vtv,, ~ (Cu.u, • C S )C S — — 1 1 C Vr 1 1 2 C ss 

which is geometrically equivalent to the following more simpler flow (by adding 
a tangential component): 

C t = C v , v , -\\\C v ,fC ss (5.10) 

Note that the flow (5.10) consists of two orthogonal diffusion terms. The 
first term C VfVf is stable as it represents a forward diffusion along the homotopy, 
while the second term — \\C Vm || 2 C SS is an unstable backward diffusion term along 
each curve. Indeed, numerical experiments show a behaviour that parallels the 
phenomenon described in §4.3.1. 

5.2 Conformal Versions of H° 

We now define the conformal H9 energy (when the conformal factor is a 
function of the arclcngth L of each curve) as 

E rj) (t) = [ <f>{L) [ \\C v ,\\ 2 dsdv (5.11) 
Jo Jo 

Once again we compute (in the appendix) the derivative of this energy along a 
time varying family of homotopies C(u, v, t). 

E'^t) = - J a J o Cf{2^'L v ,C v , + 2<f>C v , v „ - 2<f>{C v , v „ ■ C S )C S 
-2<j){C v , ■ C SS )C V , + (0771 + <j)'M)C ss ^ dsdv 

where 

m = \\C v ,\\ 2 and M= mds= / \\C V J\ 2 ds. 

Jo Jo 

As before, we now consider the planar case in which C Vt and C ss are linearly 
dependent and therefore {G Vt ■ C ss )C Vt = mC ss , yielding 

E'(t) = -2 j J C t -(0(a.t ) -(a.«;a)C' a )+0'i^.a.+^(^ , M-^n)C„)dfldt; 

(5.12) 

from which we obtain the following minimizing flow 

C t = 4>{C V , V , - {C VtV , ■ C S )C S ) + 4>'L V ,C V , + -((/>' M - <i>m)C ss 
which is geometrically equivalent (by adding a tangential term) to 

C t = 4>C V , V , + 4>'L v ,C Vt + if/M - <bm)C ss (5.13) 



5.3 Numerical results 
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5.2.1 Stable conformal factor 

To stabilize the flow in last equation, we look for a <fi such that 

(j>'M - (f>m > for all (s, w*) (5.14) 

or (assuming M ^ 0) 

j = (\og^Y>^ for all («,«») (5.15) 
One way to satisfy this is to choose 

771 

(log</>)' = max — = A (5.16) 

giving us 

<p = e XL (5.17) 
yielding the following flow of homotopies 

C t = e Ai (2C V , V , + 2XL V ,C V , + {XM - m)C ss ) (5.18) 

The choice of having <f> = e XL satisfies (5.5.*); and it agrees also with the 
discussion in §5.4, that hints that the energy E(C) associated to the metric 
may be lower semi continuous when </>(c) > len(c). The above conformal metric 
does not entail an unique Riemannian Metric on the space of curves: indeed the 
choice of A depends on the homotopy itself. 

In the numerical experiments shown in this paper, we chose A to satisfy 
(5.16) at time t = and found that this was enough to stabilize the flow up to 
convergence. However, we have no mathematical proof of this phenomenon. 



5.3 Numerical results 

We note that the minimizing flow (5.13) consists of two stable diffusion terms 
and a transport term. As such, we have the option to utilize level set methods 
in the implementation of (5.13). 

We represent the evolving homotopy C(u, v, t) as an evolving surface S(u, v, t) 

S(u,v,t) = (C(u,v),v,t) 

We then perform a Level Set Embedding of this surface into a 4D scalar function 
ip such that 

ip(C(u,v,t),v,t) = 0. 

The goal is now to determine an evolution for tp which yields the evoluton 
(5.13) for the level sets of each of its 2D cross-sections. Differentiating 

^(i)(x(u,v,t),y(u,v,t),v,t) =o) — ► Vt + VV-C t = 

where Vip = (ip x , t/) y ) denotes the 2D spatial gradient of each 2D cross-section of 
tp, and substituting (5.13), noting that N = 'Vtp/\\ V^|| , yields the corresponding 
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Level Set Evolution. 

4; =** - ij^p ( vv. ■ m + ij^jj ( v Vv*) - v« 

Note that for simplicity we have dropped the factor e XL from (5.13) since we 
are guaranteed that this factor is always positive. As a result, we do not change 
the steady-state of the flow by omitting this factor. 

If we numerically compute the geodesic of two curves Cq,Cx in figure 6, we 
obtain the geodesic, which is represented, by slicing it in figure 7 on the next 
page and as a surface in figure 8 on the facing page. 



Curve Co Curve c\ 

Figure 6: Curves Co and c\ 



5.4 Example 

This example shows why we think that the conformal energy E(C) may be lower 
semicontinuous (on planar curves) in the case when <f>(c) > lcn(c). 

Fix < e < 1/2 and A > 0. 

Suppose 

{(u, uX) if u G [0, e] 

(u,(2e-u)X) ifue[e,2e] (5.19) 
(u, 0) if u £ [2e, 1] 

is the curve in figure 9 on page 46. 

We define the homotopy C : [0, 1] x [0, 1] -> R 2 by 

C(u,v) = c(u) + (0,v) 

Let with C(u,v) = (u, v) be the identity. 

We may tessclate, as explained in point 3 in 4.3.1, to build a sequence of 
homotopies Cu such that Ch — *h C m L°° and 

d u ,vC h -±* h d UjV C weakly* in L°°, 




Figure 8: Surface of homotopy 
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Figure 9: The curve from eq. (5.19) 



and E{C h ) = E{C). 
Now we compute 



len(c) = 1 + 2e(Vl + A 2 - 1) = 1 + 2ea 



where we define a = \/l + A 2 — 1 for convenience. Note that a > 0. 
We compute the energy (using the identity (4.10)) 

E(C h ) = E{C) = [ [ f^-dudv = 



o Jo \duC\ 

(1 - 2e)0(C) + 2e-^£L = 0(C) (l-2e + 2s- 1 



VTTa 2 V " + 1 

^ (c) ( 1 _ 2e _^_ l > 

a \ „ / a 2ea 2 

> (1 + 2ea) 1 - 2e = 1 + 2e [a 



a + 1 J \ a + 1 a + 1 

1 + 2 £ a2 - W = 1 + 2 £a 2 I^ > 1 = J^(C) 
a + 1 a + 1 ~ v ; 



A More on Finsler metrics 



We now provide two more results on Finsler metrics 1.2, for convenience of the 
reader. The first result explains the relationship between the length functional 
len(£) and the energy functional E(£) 

Proposition A.l Fix x, y in the following, and let A be the class of all locally 
Lipschitz paths 7 : [0, 1] — > M connecting x to y. 

These are known properties of the the length and the energy. 

• If 7 are locally Lipschitz and is a monotone continuous function such 
that £ = 7 o and 0(0) = 0, 0(1) = 1 then Len7 = Lcn£. 

• In general, by Holder ineguality, £^(7) > Len(7) 2 . 

• If 7 provides a minimum of min^i?(7) ; then it is also a minimum of 
min_4Len(7) in the same class, E{^) = Len(7) 2 ; and moreover |7(i)| 7 (t) 
is constant in t, that is, 7 has constant velocity. 
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• If £ provides a minimum of min_4 Len(7), then there exists a monotone 
continuous junction <f> and a path 7 such that £ = 70^ and 7 is a 
minimum of mm a Ei^j). ( 48 ) 

Proof. By 2.27, 2.42, 3.8, 3.9 in [Men] and 4.2.1 in [ATOO]. □ 



This second result is the version of the Hopf-Rinow theorem 1.6 to the case 
of generic metric spaces. 

Theorem A. 2 (Hopf-Rinow) Suppose that the metric space (M,d) is locally 
compact and path-metric, then these are equivalent: 

• the metric space (M, d) is complete, 

• closed bounded sets are compact 

and both imply that any two points can be connected by a minimal length geodesic. 

A proof is in §1.11 and §1.12 in Gromov's [Gro99], or in Theorem 1.2 in [Men] 
(which holds also in the asymmetric case). 



B Proofs of §2.1.4 

Let Z be the set of all 6 £ L 2 ([0,2ir]) such that 6(s) = a + k(s)ir where k(s) £ Z 
is measurable, (a = 2ir — J fc), and 

\{k(s) = mod 2} I = \{k(s) = 1 mod 2} I = n 

Z is closed (by thm. 4.9 in [Bre86]). We see that Z contains the (representations 
6 of) flat curves £, that is, curves £ whose image is contained in a line; one such 
curve is 



&(*) = 6(«) 



3/y/2 S £ [0, 7r] 6 — \ 71 ^ S ^ ^' ^ 

(2tt-s)/V2 se(7r,27r]' 1 3tt/2 sG( 7T, 2tt] 



We provide here the proof of 2.12: M \ Z is a manifold. 

Proof. Indeed, suppose by contradiction that V^>i, V^a, V</>3 are linear depen- 
dant at 8 £ M , that is, there exists a £ R 3 , a 7^ s.t. 

ai cos(0(s)) + 02 sin(0(s)) + 03 = 

for almost all s; then, by integrating, 0,3 = 0, therefore a\ cos(6'(s))+a2 sin(^(s)) = 
that means that 6 £ Z. See also §3.1 in [KSMJ03]. □ 

This is the proof of 2.15: 



( 48 )in writing £ = 7 o </>, we, in a sense, define 7 by a pullback of £: see 2.29 in [Men]. Note 
that we could not write, in general, 7 = £ 4>~ 1 - indeed, it is possible that a minimum of 
min.4 Len(7) may stay still for an interval of time; that is, we must allow for the case when <j> 
is not invcrtiblc. 
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Proof. Fix O e M\Z. Let T = T eo M be the tangent at 6 . T is the vector 
space orthogonal to V</>j(#o) for i = 1, 2, 3. Let = ei(s) £ L 2 fl be near 
0i(#o) ha so that the map (x, y) : T x R 3 — * L 2 

3 

(x, y) h-> 9 = 6»o + a: + ^ (S.l) 
i=i 

is an isomorphism. Let M ' be M in these coordinates; by the Implicit Function 
Theorem (5.9 in [Lan99]), there exists an open set U' C T, € f/', an open 
V' C R 3 , € V, and a smooth function /:[/—> Ft 3 such that the local part 
M' n ([/' x V) of the manifold M' is the graph of y = f(x). 

We immediatly define a smooth projection ir : U' x V' — * M' by setting 
ir'(x,y) = (x,f(x)); this may be expressed in L 2 ; let (x(9),y(6)) be the inverse 
of (B.l) and [/ = (£/'); we define the projection tt : U — > M by setting 

3 

Then 

3 

7r(0)(s)-0(s) = ^e i (s)a j , o,i=(fi(x(6)) —yi) € M (B.2) 
t=i 

so if 0(s) is smooth, then ir(6)(s) is smooth. 

Let 9 n be smooth functions such that 6 n — > in L 2 , then 7r(0„) — > $o; if 
we choose them to satisfy n (27r) — 6* n (0) = 27r/i, then, by the formula (B.2), 
-k{9)(2tt) — 7r(0)(O) = 2irh so that 7r(0„) S M and it represents a smooth curve 
with the assigned rotation index h. □ 



C E N is ill-posed 

The result C.l following below was inspired from a description of a similar 
phenomenon, found on page 16 of the slides [Mum] of D. Mumford: it is possible 
to connect the two segments cq(u) = (u, 0) and c\(u) = (u, 1) with a family of 
homotopies C k : [0, 1] x [0, 1] -> IR 2 such that E N (C k ) ->k 0. We represent the 
idea in figure 10 on the facing page. 

We use the above idea to show that the distance induced by E N is zero. ( 49 -* 

Proposition C.l (E N is ill-posed) Fix cq and c\ to be two regular curves. 
We want to show that, Ve > 0, there is a homotopy C connecting cq to c\ such 
that E N {C) < e. 

1. To start, suppose that c\ is contained in the surface of a sphere, that is, 
|ci| = 1 is constant. Suppose also that |ci| = 1. Consider the linear 
interpolant, from the origin to C\ : 

C{6,v)=vcx{6) 

The image of this homotopy is a cone. 

( 49 )\Ve have recently discovered an identical proposition in [MM]: we anyway propose this 
proof, since it is more detailed. 
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\ir N d v C\ 2 \d e C\ = 
= \d v C\ 2 \d e C\sm 2 (a) ~ 
\d v C\ 2 
~ |9eC| 

case u £ [1/2, 1] 

Figure 10: Artistic rendition of the homotopy C k , from [Mum] 



We want to play a bad trick to the linear interpolant: we define a homotopy 
whose image is the cone, but that moves points with different speeds and 
times. Let e = ir/k in the following; we define the sawtooth Z : S 1 — * [0, e] 

if0e[Q,e] 

(2e-0) if6e[e,2e] 

Z{6) = {(9-2e) ifOe [2e,3e] 

(Ac -6) j/0€[3e,4e] 



(note that Z{0) + Z{0 + e) = e, Z(0) = Z(-0) ) 
Let 

2v 

C k {6,v) = c x {0)—Z(6) 
e 

forve [0,1/2], and 

C fc (M) = ci(0) (l-2i^Z(0 + e) 
for ve [1/2,1]. 

TVie energy E(C k ) is splitted in two parts for v and in 2k equal parts in 9, 
so we compute the energy only for two regions, and then multiply by 2k. 



In region 6 £ [0, e] v <G [0, 1/2], we have 



C k (d,v) = Cl (9)2-v6 
e 



d v C k =c x {6)2U 

d e C k =ci (0)2-0 + Cl (0)2- 
e e 

and 

|9 C fe | 2 =4« 2 i(0 2 + l) 



C E N IS ILL-POSED 



Since 



then 



\tnv\ 2 = \v - (v, T)T\ 2 = \v\ 2 - ((v, T)f 



\-K N d v C k \ 



e 



l0 2 (| Cl (0)| 2 -((T, Cl (0))) ; 



-e 2 (i~(d e c, Cl (0)) 2 



\d B C\< 



i-( Cl {e)2-, Cl {e) 



2 1 



\d 9 C\< 



1 - 4^rU 



1 , 1 



e- \d e C\< 



1 - 



1 



4 M 1 



i + e 2 J e 2 i + e 2 

so that the energy for the part v G [0, 1/2] of the homotopy is 

,1/2 ,e 

Jo Jo 

\d e C k \ d6 dv = 

JO JO fc ~ J--I-CT" 

1 r 1 ^ r * . ! / — 

16k- 
= 32k 



Jo Jo e 2 l + < 

/•1/2 /-e , / ; 

4w 2 — (6» 2 + 1) d6 dv 



JO 

i r 1/2 r<* i 



Jo 



1 + 6 2 



v^J(9 2 + l) dO dv 



i i f £ i 
= 32fc^- / ft A <W< 
e 3 8 io VWTTj ~ 

similarly in region 9 G [0, e] w € [1/2, 1] we /lave 



l-2l— ^(e- 



we implicitely change variable 9 t— > 9 — e, v t— > v — 1 £0 wrzie 
C fc (0,u) = ci(0) (i- 2 ^) 
ff/ws means that we will integrate on 9 G [—£,0] i> G [—1/2,0],). TTien 

1 



0„C fc = - Cl {9)29- 



and 



d e C k =c 1 {9) (l-2-^9)- Cl {9)2 V - 
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and 

\d g C k \ 2 = (l - 2 V) 2 + 4^ = I ((e - 2^) 2 + 4« 2 ) 



|7Tjv5„Cfe| 



= -9 2 {\c 1 {9)\ 2 ~{T,c 1 {9)Y 



-9 2 l-( Cl (0)2-,ci(e) 



|9 e C fe | 



1 - 4- 



1_ 

V |9eC fc | 



(e - 2w6 



(e - 2v6») 2 + 4w 2 



./Vote £/ia£ 



2c J 



(e - 2v9Y + Av 2 > e 2 (l + 2vf + Av 2 > ^ (C.l.*) 

(the positive minimum is reached at v = — e 2 / (2 + 2e 2 ), 9 = — e) 
Since 



\7T N d v C k \ 2 \d C k \ = ^e 2 



(e - 2v6 



(e - 2u(9) 2 + 4v 2 



-^( (e- 2u(9) 2 + Av 2 



(e - 2t;£ 



£ 2 (1-2wt) 2 



v /(e-2v9) 2 +4v 2 



y/e 2 (l-2vr) 2 + Av 2 



where t = 9/e so the energy for the part v € [1/2, 1] 0/ £/ie homotopy 
becomes 



-1/2 J-e 



2fc 



1/2 







e 2 (l-2wr) 5 



8tt 



4t 2 

1 v / e 2 (l-2wr) 2 +4w 2 
(l-2vr) 2 



1/2 



V(l - 2wr) 2 + 4u 2 /e 



dr dv = 
dr dv 



since by (C.l.*) the integrand is continuous and positive, and it de- 
creases pointwise when e — -> ; then so E N (Ck) — ► as e — > (by 
Beppo-Levi lemma, or Lebesgue theorem). 



2. as a second step, consider a generic smooth curve c\; we could approximate 
it with a piecewise smooth curve c[ , where each piece in C is either 
contained in a sphere, or in a radious exiting from 0: by using the above 
homotopy on each spherical piece, and translating and scaling each radial 
piece to 0, we can have a homotopy 



52 



D DERIVATION OF FLOWS 



3. then, given two generic smooth curves cq,c\, we can approximate them as 
above to obtain c' a ,c[, and build a homotopy 

c — ► c — ► — * c[ — > ci 

this final homotopy can be built with small energy 

Remark C.2 More in general, if a > [3 > , then the energy 

J \n N d v C\ a \Cf 
is ill-defined, as is shown in 4.. 3.1. 

Proposition C.3 (E is ill-posed) Consider the energy E(C) associated to 
the metric (2.6). Fix Co and c\ to be two regular curves. Ve > 0, there is a 
homotopy C G C connecting Co to c\ such that E{C) < e. 

Proof. Consider a homotopy defined as in the previous proposition; if we mollify 
it, we can obtain a regular homotopy C such that \\C — C||^i,3 and \\C — C\\oo 
are arbitrarily small; then we can reconnect C to Co to Ci with a small cost, to 
create C": by some direct computation, E{C") < 2e. 

By using prop. 3.8 on C" we obtain a C such that E N (C) = E N (C"), and 
since tt^C = 0, 

E{C) = E N (C) = E N (C") < 2e 

□ 

D Derivation of Flows 

In this section we show the details of the calculations of the minimizing flows 
for both the H° and conformal energies. 

D.l Some preliminary calculus 

First we develop in the following subsections some of the calculus that we will 
need to work with the geometric parameters s and introduced in section 5. 

D.l.l Commutation of derivatives 

Note that the parameters s and v* do not form true coordinates and therefore 
have a non-trivial commutator. The third parameter t will come into play later 
when we consider a time varying family of homotopies C(it, v,t) and take the 
resulting time derivative of either the H° or the conformal energy along this 
family. 



D.l Some preliminary calculus 
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d d d f d C u -C v d 

dt dv* dt \dv C u ■ C u du 

d_d_ C u -C v d d 

dt dv C u ■ C u dt du 

d C u -C v d\ d_ 
dv C u ■ C u du J dt 



C U t ' C V + C U ' Cy 



(C u ■ C v )(C u t ■ C u ) \ d 



C u • C u 



C u ■ (Cf. 



d d d 
(C s ■ Ctv, + Cts • Cu,) 



dv* dt 



ds 



C u -C v n 



(c u ■ c u y 

r (n - Cu-Cy f 

^tu \ 







C u • C u 



du 
(D.l) 



d_d_ d ( 1 d \ 1 d d C ut ■ C u d 

di~d~s Wt\\\C u \\lf^)~\\Cu\\~dt~du~~ \\C U \\ 3 du~ 



os at os 



(D.2) 



A A 9 ( 1 d \ c c d ( 1 d \ 

dv* ds dv \||C U || du) " ds \I|C U || du J 

i d d CuvCu d 

dv du ||C„|j 3 du v ' s ds ds 
r ^__r rll. 

OS OV OS OS OS 

~ tt [ t: C v ■ C s — J + C v ■ C ss — 

(Is \ov Os J os 



d d d 

t: — ^ 1 - C Vt ■ C ss — 

os ov* Os 



(D.3) 



D.l. 2 Some identities 

Here we write down some useful identities regarding various derivatives of the 
homotopy with respect to the geometric parameters s and v* . 



1. C s ■ C s = 1 

2 . C s ' Cy = 

^* Cy^s ' C S — ~C SS ' Cy^ 

5- C sv ^ • C s = 

6. C ss ■ C s = 



Cv*v* ' C s 



D.l. 3 Commutation of derivatives with integrals 

Finally, we write down how to commute derivatives and integrals when differ- 
entiating with respect to t or v*. 
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dt 



fds = 



d_ r 1 

diJo 



f\\C u \\du= / f t \\C u \\+f(C U fC s )du= / f t + f(C ts -C s )ds 
Jo Jo 

ft- fs{CfC s )- f(CfC aa )ds (D.4) 



fds 



d_ r 1 

9v Jo 

L 

fv — fs(C v ■ C s ) — f(C v ■ C ss ) ds 



f\\C u \\du= I f v \\C u \\+f(C uv -C s )du = J^ f v + f(C vs -C s )ds 

fv,-f(C v ,-C ss )ds (D.5) 



D.1.4 Intermediate Expressions 

The last step, before begining the flow calculation will be to introduce a few "in- 
termediate" expressions that will help keep the the expressions in the upcoming 
derivations from becoming too lengthy. 



m 



m 






(D.6) 


m t 


— 2 C Vr t ■ C Vt = 


2 Ct Vt , ■ Cy t 


(D.7) 


m s 


2 C VvS ' 


2 Ogy^ ' Cy^ 2 Cy^y^ ' C S 


(D.8) 




— 2 C v ^ v# ■ C v ^ 




(D.9) 




2 Cy^y^y^ ' Cy^ 


~t~ 2 Cy^y^ ' Cy^y^ 


(D.10) 



M v . v . 



M = [ Cy,-C v ,ds (D.ll) 
Jo 

M t = [ 2C tv , -C v , +2{C V , V , -C s )(C t ■C s )-m(C t -C ss )ds (D.12) 
Jo 

My, = I 2 C v , v , -Cy,-m (Cy, ■ C ss ) ds (D.13) 



/ 2C( (ttt> • C Vt + 2 C v t „ t • C VtVt 771 (C^ ' C SSVti ) (D.14) 
Jo 

— m (C 

^ C*v#v*v* ' "I" ^ ■ C u>kVi)i ~\~ 2 (Csi!^ ' Cy* ) -|- 777, (C S v* ' Csv* ) 

~tti{C v ^v^ 'Css) 4 (C v#Vl> ■ C V ^(C V ^ • C ss ) ds 
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L = [ ds (D.15) 
Jo 



L 



-- I -C t -C„da (D.16) 
Jo 

L v , = I -C v ,-C„d8 (D.17) 



/ C'u*u* ' C ss Cv, ' Cssv t ~\~ (Cu„ ■ Css) ds (D.18) 
Jo 



L 

C sv# ' Csv* ds @v*v* • C s 







D.2 H° flow calculation 



We are now ready to begin the flow calculation. We'll start with the case of 
the H° energy in this subsection and then proceed to the conformal case in the 
followin subsection. 



We begin by considering a time- varying family of homotopics C(u,v,t) : 
[0, 1] x [0, 1] x (0, oo) — > E" and write the H° energy as 



E(t) = f f 
Jo Jo 



|C„„|| dsdv= I Mdv 
'o 



(D.19) 
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Then the variation of E is 

r 1 pi pL 

E'{t) = M t dv= / 2 C tv , -C v ,+2 (C VmVm ■ C s )(C t ■ C.) - m (C t ■ C ss ) ds dv 
Jo Jo Jo 

= J J 2 (c tv - {C v ■ C a )C ts ) ■ C v , ds dv + J J C t -{2 (C v , v , ■ C S )C S - m C ss ) ds d« 

= 2 f f C tv -C v ,\\C u \\dudv 
Jo Jo 

+2 f f (C v ■ C s ){C t ■ C v ,„) + {C vs -C s + C v - C ss )(C t ■ C v , ) ds dv 
Jo Jo 

+ J J C t -{2 (C„.„. ■ C S )C S ~ m C ss ^j ds dv 

= 2 / f -{C t - C V ,)(C UV ■ C s ) - C t ■ C V , V \\C U \\ dudv 
Jo Jo 

+2 / f (C v ■ C s ){C t ■ C v ,„) + {C vs -C s + C v - C ss )(C t ■ C v , ) ds dv 
Jo Jo 

+ J J C t -(2 {C v , v , ■ C S )C S - m C ss ) ds dv 
= 2 J J C t • ( - C v , v + (C v ■ C S )C V , S + (C u • C„)C„.) ds dv 

+ fo Jo Ct '( 2 i ' Cv ' V ' ' ° S ' )Cs " m CsS ^ ^ ^ 
= - j Cf (2C V , V , - 2(C V , V , ■ C S )C S - 2(C„. ■ C„)C V , ■ ,n(\. ) dsdv 

In the planar CclSG, Cy cLIld Css 

are linearly dependent (as both are orthog- 
onal to C s ) which means that 

(C Vt ■ C ss )C Vt = (C Vt ■ C Vt )C ss = TnC ss (-D.20) 

and therefore 

E\t) = ~ J J o Cf (2(C V , V , - {C v , v , ■ C S )C S ) - mC ss ) dsdv (D.21) 

by which we derive the minimization flow 

Ct = 2{C VtVt {C Vt , v * ' Cs)C s ) mC ss 

D.3 Conformal flow calculation 

We now define the energy as 

p 1 t* L p 1 

E 4> {t) = I <f){L) \ \\C V ,\\ 2 ds dv = / <f)Mdv (D.22) 
Jo Jo Jo 
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and compute its derivative as 



E'(t) = I <f>'L t M + cj)M t dv 



o 



jf (-<t>'M jf C t -C ss ds + <t> 2 C to , • (7 fl . + 2 (C.„. • C s )(C t ■C a )-m (C t ■ C ss ) ds^j dv 



1 rL 

2<b(C tv - (C v ■ C s )C ts j ■ C v . dsdv 



O JO 



+ J o J o Ct - (2<t>(C v , v , ■ C S )C S - (0m + 0'M)C ss ) dsdv 

= 2 f [ <f>C tv -C v ,\\C u \\dudv 
Jo Jo 

+20 / / (C v ■ C s ){C t ■ C v , a ) + (C vs -C s + C v - C ss )(C t ■ C v , ) ds dv 



o jo 

C t ■ [2<j>{C v , v , ■ C S )C S - (<j>m + <j>'M)C ss ) dsdv 
o Jo v ' 

2 f f -((f>Cf C V ,){C UV ■ C s ) - C t ■ (4'L V C V , + (f>C v , v )\\C u \\ dudv 
Jo Jo 

+2<j> [ [ (C v ■ C s )(C t ■ C v , s ) + {C vs -C s + C v - C ss )(C t ■ C v , ) ds dv 
Jo Jo 

1 rL 

C t ■ [2(j)(C v , v , ■ C S )C S - (<fm + <j>'M)C aa \ dsdv 



10 Jo 

= 2 



to Jo 

fl rL 



C t ■ ( - <f/L v C v , -^C v , v + <t>{C v ■ C S )C V , S + 4>{C V ■ C„)C„.) dsdv 
C t ■ (2<j){C v , v , ■ C S )C S - (0m + <j)'M)C s .}j dsdv 
C t ■ h(j)'L C Vr + 2(j>C VtVt - 2(j>(C VtVt ■ C S )C S - 2<j){C v , ■ C ss )C Vt + (<f>m + (j>'M)C ss ) 



In the planar case, C Vt and C ss are linearly dependent (as both are orthog- 
onal to C s ) which means that 

{C v , ■ C ss )C Vt — [C Vt ■ C Vr )C ss = mC sa (-D.23) 

and therefore 

E'(t) = - jf jf C t -{2<f>{C v „ v ,-{C v , v ,-C a )C a )+2<l>'L v ,C v M$ 

{DM) 



to Jo 

which entails the flow 



C t = 2$(C V , V , - (C v „ v , ■ C S )C S ) + 2<j>'L v ,C v , + [4>'M - 4>m)C a 



58 CONTENTS 

Contents 

1 Introduction 1 

1.1 Riemannian and Finsler geometries 3 

1.1.1 Riemannian geometry 4 

1.1.2 Geodesies and the exponential map 4 

1.1.3 Submanifolds 5 

1.2 Geometries of curves 6 

1.2.1 Finsler geometry of curves 7 

1.3 Abstract approach 7 

2 Examples, and different approaches and results 10 

2.1 Riemannian geometries of curves 10 

2.1.1 Parametric (non-geometric) form of H° 10 

2.1.2 Geometric form of H° 11 

2.1.3 Michor-Mumford 12 

2.1.4 Srivastava et al 13 

2.1.5 Higher order Riemannian geometry 14 

2.2 Finsler geometries of curves 15 

2.2.1 L°° and Hausdorff metric 15 

2.2.2 L 1 and Plateau problem 16 

3 Basics 16 

3.1 Relaxation of functionals 16 

3.2 Curves and notations 17 

3.2.1 Homotopies 19 

3.3 Preliminary results 20 

3.4 Homotopy classes 23 

3.4.1 Class C 23 

3.4.2 Class F of prescribed-parameter curves 23 

3.4.3 Factoring out reparameterizations 24 

4 Analysis of E N (C) 24 

4.1 Knowledge base 25 

4.2 Compactness 26 

4.3 Lower semicontinuity 27 

4.3.1 Example 30 

4.4 Existence of minimal geodesies 33 

4.4.1 Michor-Mumford 35 

4.5 Space of Curves 37 

4.6 The pulley 37 

5 Conformal metrics 38 

5.1 The Unstable H° Flow 41 

5.1.1 Geometric parameters s and u* 41 

5.1.2 H° Minimizing Flow 41 

5.2 Conformal Versions of H° 42 

5.2.1 Stable conformal factor 43 

5.3 Numerical results 43 

5.4 Example 44 



LIST OF FIGURES 59 

A More on Finsler metrics 46 

B Proofs of §2.1.4 47 

C E N is ill-posed 48 

D Derivation of Flows 52 

D.l Some preliminary calculus 52 

D.l.l Commutation of derivatives 52 

D.l. 2 Some identities 53 

D.l. 3 Commutation of derivatives with integrals 53 

D.l. 4 Intermediate Expressions 54 

D.2 H a flow calculation 55 

D.3 Conformal flow calculation 56 

List of Figures 

1 Rcparamctcrization to irfd v C = 22 

2 Stretching the parameterization 22 

3 Tesselation of homotopy C to form Ch 31 

4 The pulley Ch, in case h = 5 38 

5 Intersections induced by flow 39 

6 Curves Co and ci 44 

7 Slices of homotopy 45 

8 Surface of homotopy 45 

9 The curve from eq. (5.19) 46 

10 Artistic rendition of the homotopy Ck, from [Mum] 49 



References 

[AFP00] Luigi Ambrosio, Nicola Fusco, and Diego Pallara, Functions of 
bounded variation and free discontinuity problems, Oxford Mathe- 
matical Monographs, The Clarendon Press Oxford University Press, 
New York, 2000. MR MR1857292 (2003a:49002) 

[AT00] L. Ambrosio and P. Tilli, Selected topics in "analysis in metric 
spaces", appunti, edizioni Scuola Normale Superiore, Pisa, 2000. 

[Atk75] C. J. Atkin, The Hopf-Rinow theorem is false in infinite dimensions, 
Bull. London Math. Soc. 7 (1975), no. 3, 261-266. MR MR0400283 
(53 #4118) 

[BCS] D. Bao, S. S. Chern, and Z. Shcn, An introduction to Riemann- 
Finsler geometry, (October 1, 1999 version). 



[Bre86] 



H. Brczis, Analisi funzionale, Liguori Editore, Napoli, 1986, (italian 
translation of Analyse fonctionalle, Masson, 1983, Paris). 



GO 



REFERENCES 



[But89] G. Buttazzo, Semicontinuity, relaxation and integral representation 
in the calculus of variation, scientific & technical, no. 207, Longman, 
1989. 

[CFK03] G. Charpiat, O. Faugeras, and R. Kerivcn, Approximations of shape 
metrics and application to shape warping and empirical shape statis- 
tics, INRIA report 4820, 2003. 

[Dac82] Bernard Dacorogna, Weak continuity and weak lower semicontinuity 
of non-linear functionals, Lecture Notes in Mathematics, vol. 922, 
Springer- Verlag, 1982. 

[EE70] J. Eells and K. D. Elworthy, Open embeddings of certain Banach 
manifolds, Ann. of Math. (2) 91 (1970), 465-485. MR MR0263120 
(41 #7725) 

[Eke78] Ivar Ekeland, The Hopf-Rinow theorem in infinite dimension, J. Dif- 
ferential Geom 13 (1978), no. 2, 287-301. 

[Fom90] A.T: Fomcnko, The Plateau problem, Studies in the development of 
modern mathematics, Gordon and Breach, 1990. 

[Gro99] M. Gromov, Metric structures for Riemannian and non-Riemannian 
spaces, Birkhauser, 1999. 

[Kli82] Wilhelm Klingenberg, Riemannian Geometry, W. de Gruyter, 
Berlin, 1982. 

[KSMJ03] Eric Klassen, Anuj Srivastava, Washington Mio, and Shantanu Joshi, 
Analysis of planar shapes using geodesic paths on shape spaces, 2003. 

[Lan99] Serge Lang, Fundamentals of differential geometry, Springer, 1999. 

[Men] A. C. G. Mennucci, On asymmetric distances, preprint, 
http : //cvgmt . sns . it/papers/ and.04/. 

[MM] Peter W. Michor and David Mumford, Rie- 
mannian geometris of space of plane curves, 
http://front.math.ucdavis.edu/math.DG/0312384. 

[Mum] D. Mumford, Slides of the gibbs lectures, 

http://www.dam.brown.edu/people/mumford/Papcrs/Gibbs.pdf. 

[Sha49] C. E. Shannon, A mathematical theory of communication, The Math- 
ematical Theory of Communication, University of Illinois Press, 
1949. 

[Sim83] L. Simon, Lectures on Geometric Measure Theory, Proc. Center for 
Mathematical Analysis, vol. 3, Australian National University, Can- 
berra, 1983. 

[YM04] A. Yezzi and A. Mennucci, Conformal riemannian 
metrics in space of curves, 2004, EUSIPCO04, MIA 
http://www.ceremade.dauphine.fr/ cohen/mia2004/. 

[You98] Laurent Younes, Computable elastic distances between shapes, SIAM 
Journal of Applied Mathematics 58 (1998), 565-586. 



